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Te role of James Clerk Maxwell as the father of radio science is well known and his work has been 
honoured by this Institution with the establishment in 1951 of the Clerk Maxwell Memorial Lectures. Since 
1945 a senior premium awarded by the Institution has also borne his name. 

Maxwell was born on 13th June 1831 in Edinburgh and this year therefore marks the 150th anniversary of 
his birth. His early life was spent at Glenlair, Kirkcudbrightshire, and it was there he was buried after his death 
in Cambridge on 5th November 1879 at the age of forty-eight years. He entered Edinburgh University in 1847 
and after three years went on to Cambridge, becoming a student at Trinity College and passing the 
mathematical tripos with distinction in 1854. 

For the next two years Maxwell was at Cambridge, lecturing in hydrostatics and optics in which he carried 
out many experiments, and ‘working away at electricity’. In 1856 he was appointed to the Chair of Natural 
Philosophy at Marischal College, Aberdeen where he remained until 1868 when he was appointed to the Chair 
of Natural Philosophy at King's College, London. During the next five years he worked mainly on colour but 
also laid the foundations for his later, great work, the ‘Treatise on Electricity and Magnetism’. From 1865 he 
was in retirement at Glenlair and apart from occasional special lectures and acting as examiner for the 
mathematical tripos at Cambridge, his main energies were devoted to his ‘Treatise’ which was published in 
1873. 

Physics at Cambridge was on the threshold of great changes at this time, and in 1870 the Duke of 
Devonshire, then Chancellor of the University, offered to build and furnish a Physical Laboratory. In 1871 
Maxwell was appointed to the new Chair of Experimental Physics and designed and supervised the erection of 
the new laboratory which was formally opened in June 1874. 

For the next four years, Maxwell's most important and continuing work was superintending experimental 
work in the Cavendish Laboratory. His last public lecture, the Rede Lecture at Cambridge in 1878, was entitled 
"On the Telephone’, the invention of which he considered to be a ‘material symbol of the widely separated 
departments of human knowledge', and he spoke of the great benefits which could result from the cross- 
fertilization of the sciences. In a fine tribute on the occasion of the centenary of his birth, Max Planck said of 
him ‘It was his task to build and complete the classical theory, and in so doing he achieved greatness unequalled. 
His name stands magnificent on the portal of classical physics, and we can say this of him: by his birth James 
Clerk Maxwell belongs to Edinburgh. by his personality he belongs to Cambridge, by his work he belongs to 
the whole world.” 

F.W.S. 


The 10th Clerk Maxwell Memorial Lecture, entitled "The Complete Physicist’ , was given at the recent 
Institution Conference in Leeds by Professor R.V. Jones, C.B., F.R.S., who held the Chair of Natural 
Philosophy at the University of Aberdeen from 1946 until his retirement this year. The lecture, which 
assesses Maxwell's many contributions to science, will be published in the November/December 1981 issue 
of The Radio and Electronic Engineer. 
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Letters to the Editor 


From: B. Priestley, B.Sc., C.Eng.. M.I.E.R.E. 
A. H. Harrold (Graduate) 
G. E. Smythe (Associate) 
T. Lomas, E.R.D. Dip.E.E., 
M.LE.R.E. 


C.Eng, M.LE.E., 


Oscillator Design 


I was gratified to note from the May issue of The Radio and 
Electronic Engineer (p. 203) that others share my concern over 
the state of oscillator design, and I would like to add the 
following comments. 

First to Mr Mitson: | did not refer solely to oscillators 
produced for biological telemetry although 1 did include some 
of the horrible examples sold to unsuspecting biologists by so- 
called engineers. If this makes a biologist feel he can do no 
worse himself, it underlines the need for liaison indicated by 
Mr Mitson (plus a gentle pointer on the mis-use of the term 
engineer?). 

Secondly to Messrs Kentley and Rogers: Although the 
quartz crystal does complicate oscillator design, the basic 
principles are exactly the same for any sinusoidal oscillator and 
it is these which are glossed over. If additional guidance is 
available in industry it can only be used over this foundation. 
(Personally I have found the attitude of crystal manufacturers 
to vary widely from ‘We consider this information will not help 
you so we won't tell you', to the very friendly and helpful.) 

However, since first writing I have discovered it is now 
possible to qualify in electronics without any knowledge of 
electromagnetic field theory! (IERE Education Committee 
please note!) In such circumstances we seem to be fighting a 
losing battle suggesting that engineers should learn basic 
principles rather than be programmed to fill a slot. 

B. PRIESTLEY 
43 Raymond Road, 
Langley, Slough, 
Berkshire SL3 8LN 
ISth May 1981 


[The Education and Training Committee is well aware of the 
present trend towards modifying qualifications and has given it 
much thought on several recent occasions. Basically the 
problem is that a three-year degree course is not long enough in 
which to cover every subject that might be considered essential. 
But the Committee is not always happy about what has been 
sacrificed as a consequence in some courses.—EDITOoR.] 


Multi-disciplinary Engineering in the Process Industries 


What a pleasant breath of fresh air greeted me as | opened my 
journal in June! There on page 255 was an article from kindred 
spirits who have at last recognized the needs of the process 
industry. | refer to the brief paper 'Multi-disciplinary 
engineering in the process industries’. 

Having started my working like in the electronics industry, 
gaining HNC and endorsements on a day-release basis, I 
followed the process instrumentation route into the process 
industry. I have over the years worked for both contractors 
and end-users and now | have my own process control 
engineering design company. 
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The need for multi-discipline knowledge is paramount in this 
field of engineering, much of which can only be gained by bitter 
experience. The engineers who work in this sector, and I am 
proud to be amongst them, have to satisfy very rigorous 
requirements before being able to follow this sometimes 
frustrating but always satisfying way of life. These are briefly: 
1. A sound technical education—not necessarily a degree but 
unfortunately this now seems to be the minimum. 

A sound practical approach to engineering, for what you 

design has to work and withstand all manner of 

environmental abuse. 

An ability to communicate with a full cross-section of 

people from directors and managers to skilled and unskilled 

artisans—without this ability you will not get off base I. 

A willingness, albeit sometimes reluctantly, to go out onto 

site in all weathers at all time of the day or night—don't 

forget this is what engineering is all about. 

. Above all other criteria, the ability to understand and 
appreciate the requirements and difficulties of engineers of 
all disciplines, and a willingness to learn from them as they 
learn from you-——there’s no place for prima donnas on site, 
teamwork is the key. 


2 


Much of what I have listed above is good sound common 
sense and it is this area of our engineers training that needs a 
lot of attention. l, on reflection, consider myself fortunate that 
I studied on a part-time basis as this permitted me to work 
alongside people at shop-floor level during the most 
informative years of my apprenticeship. I feel many graduates 
today have missed out on this experience and whilst they may 
be the managers to tomorrow they may never gain the 
'common touch'—a bridge which must be strong for our 
nation's success in the future. 

So it's back to our academics to appreciate that not all 
engineering is carried out in laboratories or factories with 
clearly defined and narrow bands of development, and it's back 
to our Institutions to recognize the vast number of sound 
practical engineers, who may not be academics, but without 
their foresight and abilities there would not be the need for 
research and development as we know it today. 

A. H. HARROLD 
"Corner House’ 
Forest Edge Road, 
Ringwood, 
Hampshire BH24 3DF 
10th June 1981 


The ‘parting shot’ in the letter from Mr. A. C. Evans (June 
1981, p. 258) caused me to reflect that his remark about the 
paucity of correspondence in the journals of the learned 
Institutions is not true of other publications such as Wireless 
World. 

Is it perhaps the case that the majority of engineers can spare 
insufficient time from their immediate problems fully to 
evaluate such erudite contributions (in the issue cited all the 
papers come from Research Establishments or Universities) to 
the extent necessary to comment on them, or even, dare one 
suggest, to understand the more specialized papers? That they 
are valuable contributions to the state of the art no one doubts 
but they seem not to inflame the passions quite as much as 
suggestions that a pair of PX 4's in class A push-pull may 
out perform any of these modern transistor output stages, or 
that Essen may supplant Einstein. 

Correspondence seems to thrive on controversy and this 
often throws up interesting points that otherwise might not 
come to light. Perhaps then you could start by inviting 
suggestions for technical topics which are open to argument, 
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such as, *A study of the optimum frequency and mode for an 

open-channel (c.b.) service in the European Common Market’. 
G. E. SMYTHE 

9 Camborne Road. 

Sutton, 

Surrey SM2 6RJ 

llth June 1981 


I am writing to say how very much I sympathize with the views 
expressed in the article "Multi-disciplinary engineering in the 
process industries' in the June issue of the Journal. There is a 
clear and imperative need in any technology based business or 
industry (I would suggest not only in the process industries) for 
people with a thorough grounding in technology plus an 
appreciation of the many other important facets of business 


e.g. marketing, sales, finance, delivery, quality, etc., if the 
business is to survive and thrive in today's highly competitive 
world. This is not to say that there will not always be a need for 
specialists but rather that the very survival of business depends 
upon those who can take a comprehensive view of business 
opportunities and threats. | would, indeed, question whether 
an engineer without an understanding of business issues other 
than technology can properly be regarded as a professional 
engineer. In drawing attention to the need for engineers to be 
trained in other than a single discipline the article has done a 
good service. 

T. Lomas 
2] Clarence Road 
St. Albans 
Herts ALI 4NP 
24th June 1981 


Announcements 


Back copies of the Journal 


Because of exceptional calls for some of this year's issues of 
The Radio and Electronic Engineer, the Institution's stocks 
have become depleted and difficulty is anticipated in fulfilling 
future new subscriptions for 1981. Any member who has no 
further need for his or her copies is therefore invited to send 
them to the Publications Department at 99 Gower Street. The 
issues required are: 
February, March, April, June 1981 
and any of these will be most welcome. 

Stocks of certain earlier issues are also exhausted notably the 
following: 

1977 January/February, March, April, August, September 

1979 February 

1980 January/February, March, July, August 

These are needed to meet orders from libraries for complete 
back sets for binding. The Institution will gladly refund 
postage. 


Subscriptions for 1982 


The subscription rates to non-members and libraries for The 
Radio and Electronic Engineer for 1982 will be as follows: 


Annual Single 
Subscription issue 
UK and Ireland £38.00 £3.50 
Overseas £43.50 £4.00 
North and South America $92.00 $8.00 
Correction 


In the paper "Integrated optics, a tutorial review', published 
in the July/August 1981 issue of The Radio and Electronic 
Engineer, the following correction should be made. 

Page 404, Sect 5.4, The final sentence in the first paragraph 
should read: This is the principle governing positive 
photoresists but, by suitable choice of polymer-solvent 
systems, a negative resist can be used in which the incident light 
beam can be employed to indicate polymerization of the 
deposited monomer. 
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Standard Frequency and Time Service 
(Communication from the National Physical Laboratory) 


Relative Phase Readings in microseconds 
NPL—Station 
(Readings at 1500 UT) 


June Droitwich 
1981 MSF 60kHz GBR 16 kHz 200 kHz* 

1 14 18-1 61-7 

2 14 18.5 61-9 

3 1:4 17-1 62:1 

4 1:5 18-6 62.3 

5 14 18:1 62:5 

6 1:6 17:8 62:8 

7 1:6 18:1 63.0 

8 1:6 17.8 63.3 

9 1:6 18-3 63-6 
10 1-6 17-7 63.8 
11 1:6 16-4 64-1 
12 1:4 15:8 64-3 
13 1:2 16:1 64-5 
14 1:2 16-1 64-7 
15 13 14:1 64-9 
16 1:2. 17:1 65-1 
17 12 16:4 65-3 
18 1:2 16:4 65-6 
19 1:2 16.3 65.8 
20 12 16:3 66-0 
21 1:2 16:3 66.3 
22 0-9 16-7 66-5 
23 1-1 167 66-6 
24 1-1 16:8 66:6 
25 11 16:3 66-9 
26 1:1 16:9 67-1 
27 1-1 16-3 67.3 
28 11 16-3 67-6 
29 0-9 16:4 67.7 
30 0:9 16:2 68-0 


Notes: (a) Relative to UTC scale (UTC, p, -Station- +10 at 1500 

UT, 1st January 1977. 

(b) The convention followed is that a decrease in phase 
reading represents an increase in frequency. 

(c) 1 us represents a frequency change of 1 part in 10'' per 
day. 

* It may be assumed that the satellite stations on 200 kHz at 
Westerglen and Burghead will follow the day to day 
changes in these phase values. 
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Function Generator Review 


During recent months, three leading companies have 
introduced new function generators. Although each has a 
different declared purpose, brief notes of their features are 
interesting as giving an insight into current practice. 


A new multi-purpose waveform generator, the 2123 has been 
launched by Marconi Instruments. This is a portable and highly 
versatile instrument, suitable for use in educational 
establishments, laboratories, production and maintenance 
areas. It costs £215 (excl. V.A.T.) 

The 2123 produces sine, square, and triangular waves at 
frequencies between 0-003 Hz and 200 kHz, with ramp 
generation between 0-003 Hz and 20 kHz. Source impedance is 
60 Q on all functions, with variable output to a maximum of 
10 V p-p, calibrated into a 600 Q load. Two auxiliary outputs 
enable a triangular wave to trigger an oscilloscope sweep, and 
TTL compatible square wave to supply up to 20 standard 
loads. 


Marconi Instruments multi-purpose waveform generator 2123 


To simplify measurement of frequency response in amplifiers 
and filters, the 2123 is fitted with an internal logarithmic sweep, 
down to one thousandth of the maxium frequency displayed on 
the dial. 

The frequency can be set by seven range pushbuttons on the 
front panel, in conjunction with a multi-turn potentiometer dial 
with three-digit read-out. Alternatively, the frequency may be 
controlled by application of a direct voltage (v.c.f.) to rear 
panel sockets; this is particularly useful for users of control or 
automatic systems. 


A high level of cost-effective performance is claimed for a 
new log sweep function generator announced by Racal-Dana 
Instruments. The Exact 500SL function generator, for which 
Racal-Dana Instruments is sole UK distributor, has a 


Exact S00SL function generator 
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frequency range of 0-005-5 MHz, with many advanced high- 
performance features, and is competitively priced at less than 
£250. 

The generator provides a haversine mode, low distortion sine 
waves (0-527), highly linear triangles—greater than 99%—and 
fast rise and fall square waves of less than 40 ns. Output 
amplitudes are 20 V peak-to-peak open circuit —10 V peak-to- 
peak into 50 Q—and may be attenuated up to 40 dB using the 
variable amplitude control and the low output socket. A d.c. 
offset control allows adjustment of waveform offset up to 
+10 V. The d.c. offset control has a detented off position when 
no offset is desired. 

Frequency may be varied manually with the multiplier 
control, externally via a v.c.f. input, or internally with a built-in 
sweep. A pushbutton selects either linear or logarithmic 
frequency variation. 

A lockout mode is set on release of the button. Sine and 
triangle signals will lock at zero volts, and the square and 
haversine signals will lock at their negative peaks. An external 
or manual trigger will produce one cycle of the selected 
waveform which stops at its normal lockout level. When the 
gate pushbutton is depressed, the external or manual signal will 
cause the generator to run until the signal goes low or the 
manual pushbutton is released. 

In addition to high, low, and v.c.f. inputs, the 500SL 
features a TTL compatible sync output, a V:f (voltage 
proportional to frequency) output, and a sweep output. It is 
housed in a rugged aluminium box providing good e.m. 
screening and weighs 2-27 kg (51b). 


A new top-level function generator has been added to Philips’ 
wide range of low frequency equipment, supplied in the 
UK by Pye Unicam of Cambridge. It is the PM 5134, a 
versatile and sophisticated research-level function generator 
which covers a wide 20 MHz bandwidth and provides sine, 
square and triangular waveforms as well as positive and 
negative going pulses and d.c. 
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Philips research-level function generator—PM5134 


Normal or crystal-controlled operation is possible—the 
latter including a frequency-lock facility and giving accuracies 
of +5 parts in 10°. Single and burst modes are also provided. 
Versatile sweep facilities allow internal or external single or 
continuous sweep with independent controls for start and stop 
frequencies. 
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Other features of the generator include amplitude and 
frequency modulation, internally or externally, with a.m. 
depth adjustable from 0 to 100% and f.m. from 0 to 10%. The 
instrument produces a high 20 V peak-to-peak output (10 V 
peak-to-peak pulses) with stepped and variable attenuation. 
Output impedance is pushbutton selectable to 50 or 600 Q. 

A three-and-a-half digit le.d. read-out displays output 
frequency and open-circuit output voltage. There are also 
separate l.e.d.s indicating frequency and voltage ranges, sweep 
modes and setting errors. The 1 mHz to 20 MHz frequency 
range is covered in ten linear sub-ranges with separate coarse 
and fine controls. Attenuation up to 60 dB is incorporated in 
calibrated steps of 3, 6, 10 and 20 dB with variable attenuation 
up to 20 dB. D.c. offset is available from —5 to -- 5 V and the 
duty cycle can be varied from 10 to 90%. 


Very stable and accurate output signals can be obtained 
using the crystal-controlled mode. A crystal lock position is 
included for repetitive and other high accuracy work. This 
facility is also available for a.m. operation. 

Sweep mode operation can be started either manually or 
electronically by remote control. The sweep period is 
continuously adjustable from 5 ms to 100 s and there are both 
hold and reset facilities. As with single-cycle and burst signals, 
the start phase is fully adjustable from —90° to +90°. 

The PM 5134 is designed for a wide range of applications in 
high level education and industrial research and development. 
It is a robust instrument, its strong metal case giving full 
electrical screening and good mechanical protection. It costs 
£1100. 


A Digital Television Studio Standard 


The Technica! Committee of the European Broadcasting Union 
(EBU) has agreed to propose a set of technical parameters as a 
digital television studio standard, in the hope that these will be 
adopted throughout the world. The decision is seen as a 
watershed in broadcasting technology. 

The 1960s left the world with a diversity of colour systems 
(PAL, SECAM, NTSC) and these will be used for broadcasting 
to the public for many years to come. Existing receivers will, 
therefore, not become obsolete in the foreseeable future. 

However, in the field of television programme production, 
digital technology offers the potential for considerable 
improvements in facilities, picture quality, and many other 
advantages. An agreed worldwide digital standard in this area 
would mean less expensive equipment, because of the 
economies of scale, and benefit international programme 
exchange. It would have important financial advantages for 
users and manufacturers, lead to improvements in the technical 
quality of programmes shown in a country with a different 
television system from the originating country, and assist the 
international exchange of technical information on television 
production. 

In order to consider not only the requirements for the 625- 
line television standards used in Europe, but also those of the 
regions where a 525-line standard is used, the EBU has 
coordinated its studies with the Society of Motion Picture and 
Television Engineers in North America. 

After considerable discussion, the EBU proposes that the 
digital standard should not be based on the PAL, SECAM or 
NTSC systems, but instead on a principle technically known as 
"component coding’ with a particular set of parameters that 
should apply both in countries using 525-line and 625-line 
television standards. 

The EBU is proposing a parameter set known at 
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113:5:675:675'. By popular convention, this describes a 
component-coded digital television system with a luminance 
signal sampled at a frequency of 13-5 MHz, and the two 
colour-diffcrence signals sampled at a frequency of 6-75 MHz. 
It is hoped that these parameters will be applicable to both 525- 
line and 625-line television systems, and hence that the number 
of samples per television line will then be very similar. Because 
of the different field rates. the number of samples per picture, 
however, will be different. 

The EBU proposal will be submitted to the forthcoming 
Final Meetings of the CCIR Study Groups in September 1981 
in Geneva, Switzerland, in the hope that it will obtain 
worldwide support. 

The need for a worldwide standard in this field was the 
subject of a recommendation by all unions at the 3rd World 
Conference of Broadcasting Unions held in Tokyo in 1980 and 
already there are positive indications of support for the EBU 
proposal from several of the other unions. In particular, the 
OTI (Organizacion de la Television Iberoamericana) has 
already conveyed its support to the EBU for the proposal. This 
is particularly noteworthy because both 525-line and 625-line 
countries exist in the OTI area. Encouraging reactions have 
also come from unions in North America, NANBA (North 
American National Broadcasters Association) and in Africa, 
URTNA (Union of National Radio and Television 
Organizations of Africa). 

The EBU Technical Committee, being the consultative 
engineering body of the EBU, has submitted its views to the 
Union's Administrative Council, and the proposal was 
expected to become the subject of an official EBU policy 
statement after the meeting of the EBU Administrative 
Council, scheduled for 22nd-25th May 1981. 
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Members' Appointments 


CORPORATE MEMBERS 


Cdr D. G. Edwards, RN (Ret) (Fellow 1964. 
Member 1950) who recently retired as Senior 
Reliability Engineer with the Ministry of 
Defence, Navy Department, has taken up 
consultancy work in the reliability and quality 
assurance field. 


E. R. Friedlaender (Fellow 1944) is now 
Deputy Mayor of the London Borough of 
Brent where he has been a councillor since 
1971. After coming to the UK in 1936 he was 
a specialist on design and manufacture of 
magnetic dust core components while before 
his retirement in 1973 he was initially a patent 
consultant and later a freelance industrial 
journalist working with European technical 
magazines. 


Vice Admiral Sir Philip Watson, K.B.E., 
M.V.O. (Fellow 1967) has recently been 
appointed Chairman of Marconi Radar 
Systems, For the past three years Sir Philip 
has been the Naval and Defence Affairs 
Adviser to GEC Marconi and a director of 
the Marconi International Marine Company. 


Cdr A. C. Cowin, M.Sc., RN (Ret) (Member 
1966. Graduate 1961) has taken up a Civil 
Service appointment as an Electrical Engineer 
with the Delence Intelligence. Staff of the 
Ministry of Defence. Cdr Cowin has served 
on the Membership Committee of the 
Institution since 1980. 


D. R. Hill, M.Phil. (Member 1980) has joined 
H. R. Smith (Microwaves) of Eynsham, 
Oxfordshire. After setting up new premises 


for the Company at Kingsland, 
Herefordshire, he will be concerned with 
research. design and development of 
microwave components and subsystems, 


notably stripline and microstrip antennas. He 
was previously Principal Antenna Engineer 


with Southern Microwave Laboratories, 
Hayling Island, Hants. 
M. Hounsell (Member 1971) has been 


appointed Manager, Programme Services and 
Engineering, South, with BBC Television. 
He joined the BBC in 1963 andn worked in 
Television Networks at the Television Centre 
in London. Since 1975 Mr Hounsell has been 
Senior Television Engineer at Plymouth. 


J. W. B. Mcellwraith (Member 1956, Graduate 
1952) who was Chief Telecommunications 
Officer, Airports and Airways Aids, with the 
National Air Traffic Services at Liverpool, is 
now Director of — Telecommunications 
(Projects) based in London. 


Lt Cdr P. Nightingale, RN (Member 1972, 
Graduate 1967) has just taken up a two-year 
appointment in the Bahamas as Training 
Design Officer with the Bahamian Defence 
Force. 


B. W. Potter (Member 1969, Graduate 1964) 
has been appointed Manager. Programme 
Services and Engineering, East with the BBC. 
A member of the Corporation since 1960, he 
has held appointments with the 
Communications Department in various 
regions, the last of these being 
Communications Manager at Manchester 
since 1979, * 


G. Salter (Member 1964) who has been with 
the BBC since 1943, latterly as Head of 
Programme Services and Engineering, Wales, 
has left the Corporation and taken up a post 
as Studio Technical Systems Advisor in the 
Network Operations and Maintenance 


Department of the Independent Broadcasting 
Authority at Crawley Court. * 


W. A. M. Smith (Member 1969) has retired 
from Cable and Wireless after twenty seven 
years of service with the company. Mr Smith 
has spent the last fifteen years in Hong Kong, 
his last appointment being Deputy General 
Manager and Divisional Manager, 
Telecommunications Service and Marketing. 


W. St John White (Member 1970, Associate 
1947) who joined the Decca Navigator 
Company in 1946 as a Development Engineer 
following wartime service with the 
Telecommunication Research Establishment, 
has been appointed to the board of Racal- 
Decca. Since 1969 he has been a director of 
Decca Navigator Company responsible for 
marketing of avionic systems. 


S. R. Tumoe (Member 1973) has been 
appointed Director of Telecommunications in 
the Sierra Leone Posts & Telecommunications 
Department. He was formerly senior engineer 
and head of the Network Planning Division. 


P. R. Thorogood (Member 1971) who has 
been with Cable and Wireless at Hong Kong 
since 1980, has been transferred to the post of 
Divisional Manager, Engineering of Bahrain 
Telecommunications Company. 


NON-CORPORATE MEMBERS 


Sir Kenneth Corfield F. Eng.. F.I. Mech. E. 
(Companion 1976), who is Chairman and 
Managing Director of Standard Telephones 
and Cables, has been awarded an Honorary 
Doctorate of Science by the City University. 
Sir Kenneth's appointment as Chairman- 
designate of the new Engineering Council was 
announced on July 30th. 


* Editorial Note Incorrect details of the 
appointments held by these Members were 
printed in the June issue of the Journal. 
Apologies are extended to the Members 
concerned and to readers for any 
embarrassment or confusion this caused. 


Obituary 


The Institution has learned with regret of the 
deaths of the following members. 


Derek Vere Staynor, M.B.E. (Member 1959) 
of Boreham Wood. Hertfordshire, died on 
15th May 1981, aged 56. Born and educated 
in India, Mr Staynor served with the Royal 
Air Force as an Instructor Mechanic during 
the war and worked initially for four years as 
an Inspector in charge of instruments. with 
Orient Airways at Karachi Airport. He then 
joined Field Aircraft Services, Bovingdon 
becoming Chief Inspector. Since 1957 he had 
been with Microwave and Electronics 
Division of Elliott Brothers, London. He was 
appointed M.B.E. in the 1971 New Year 
Honours List. 


Edwin Charles Stevenson (Member 1961) of 


Wokingham. Berkshire. died recently. aged 
50. Following National Service as an Air 
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Radar Mechanic with the Royal Air Force. 
Mr Stevenson worked for some two years 
with the Ministry of Aviation Inspectorate of 
Electrical and Mechanical Engineering before 
moving in 1955 to General Radiological as a 
Senior Electronics Test Engineer. From 
1957-60 he was Assistant Development 
Manager and Project Engineer with the 
Multitone Electric Company. In 1960 he 
joined Sangamo Weston in Enfield. 


Brian John Swainson (Member 1961, 
Graduate 1959) of Southport died on 10th 
November 1980, aged 46. Following National 
Service with the Royal Air Force, Mr 
Swainson joined the Automatic Telephone 
and Electric Company, Liverpool as 
Telecommunications Engineer in 1956. He 
remained with the company, subsequently 
Plessey Telecommunications, until the time of 
his death. For a number of years Mr 
Swainson was Honorary Secretary of the 
Institution's Merseyside Section. 


Eric Victor Taylor (Member 1964, Student 
1960) of Hutton. Essex, died on 17th March, 
aged 48. Mr Taylor undertook an apprentice- 
ship with Kelvin Hughes from 1946-51 and 
subsequently worked in the Marine Radar 
Laboratory as a Development Project 
Engineer. In 1960, following part-time study 
at East Ham College of Technology. he 
gained his HNC in electrical engineering with 
endorsements. 


Major Arthur Stanley Tate, R Sigs (Ret) 
(Member 1943) of Rainham, Kent, died in 
November 1980, aged 82. Major Tate joined 
the Royal Signals in 1916 and in 1929 was 
promoted to Foreman of Signals. In 1941 he 
was granted a commission and was for some 
time an Inspecting Officer for electrical and 
mechanical equipment. After leaving the 
army, he continued with the Inspectorate until 
his retirement in 1963. 
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New Books Received 


The following books which have been received recently have been placed in the 
Institution's Library and may be borrowed by members resident in the British 


Isles. 


Design of Continuous and Digital Electronic 
Systems 


GoRboN J.A. Birp (EMI Electronics) 


McGraw-Hill, Maidenhead, 1980. 
19 x 25.5 cm. 408 pages. £20. 
CONTENTS: Transform theory. System 
analysis. Filter transfer functions. The 


transient response of filters. Feedback and 
root locus. Active filters. Phase locked loops. 
The = transform. Sampled data systems. 
Computer simulation, Digital filters. 
Intended both for practising engineer and 
student this book deals with the methods of 
solving practical problems in electronic 
engineering through the use of transform 
theory and the concept of poles and zeros. 
Covers three broad areas: Basic theory and 
filler applications; Feedback theory and 
applications; Digital theory and applications. 


Solid State Radio Engineering 
HERBERT L. Krauss, CHARLES W. BOSTIAN 
(Virginia Polytechnic Institute and State 
University) and FREDERICK H., RAAB 
(Austin. Company). Wiley. New York, 
1980, 17 x 23.5 cm. 534 pages. £12. 
CowTENTs: Radio communication systems. 
Electrical noise. Resonant circuits and 
impedance transformation. Small-signal high- 
frequency amplifiers. Sinewave oscillators. 
Phase-locked loops. Mixers. Modulation. 
Amplitude modulation receivers. F.m. and 
p.m. receivers. Television receivers. Linear 
power amplifiers. Tuned power amplifiers. 
High-efficiency power amplifiers. C.w.. f.m. 
and a.m. transmitters. Single-sideband trans- 
mitters. 

In its analysis and design of r.f. electronic 
circuits, this book reflects the developments of 
the past decade which have initiated an 
unprecedented growth in the use of analogue 
radio systems for personal and business voice 
communications. The book is intended to be 
both a reference for the working engineer and 
a textbook for senior level students in 
electrical engineering and electrical tech- 
nology. 


Analog Signal Processing and [Instrumentation 
ARIE F. ARBEL (Haifa Institute of Tech- 

nology). Cambridge University Press, 

1980. 16 x 23 cm. £27.50. 440 pages. 
CoNTENTs: Basic and feedback amplifiers. 
Transducers. Analog signal conditioning. 
Analog signal measurement in the time 
domain. Digital-to-analog converters. 
Analog-to-digital conversion techniques., 
Converters, selected topics. Data acquision 
systems. Estimation of statistics of waveforms 
and events. Signal enhancement. Computer 
aided design. 

Intended for the electronic system designer 
who utilizes readily available integrated 
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circuit modules. Most of the material grew 
out of a course taught by the author at the 
Technion-lIsrael Institute of Technology. 


Handbook of Rectifier Circuits 
G. J. ScoLes (formery English Electric Valve 
Company). Ellis Horwood. Chichester, 
1980. 15 x 23.5 cm. 235 pages. £22. 
Contents: Rectifying devices. Half-wave, 
full-wave and bridge. Three-phase rectifier 
bridges having an output frequency of 12 
times the mains frequency. Voltage 
multiplying bridges. Biphase voltage- 
multipliers. Single-phase/three-phase hybrid 
circuits. Double-input circuits. Capacitively- 
coupled, double-driven rectifier circuits. 
Partial- bridges. Half-voltage circuits. 
Hammer circuits and anvil circuits. Double- 
output rectifier circuits. Constant-current 
techniques. Voltage control by grids and 
gates. Techniques for generating ultra-high 
voltages. Multi-stage  voltage-multiplying 
circuits. The tuned-bridge rectifier and its 
application to voltage-multiplying circuits. 
Switched, tuned, voltage-multiplier circuits. 
Further possible multiplying circuits. Three- 
phase transformer connection, mainly for 
bridge rectifier circuits. Simple smoothing 
circuits. The switching of transformers. 
Circuit component rating. 

The handbook first describes the relatively 
few rectifier circuits in common use and then 
shows how these can be combined and 
adapted so as to produce a large number of 
‘brand new! circuits of wide application 
wherever rectified power supplies are 
required. 


Electronic Devices 

F. R. Connor. Edward Arnold, London. 
1980. 13.5x21.5 cm. 121 pages. £3.95. 
(Paper bound). 

Contents: Atomic theory, Semiconductor 

theory. Solid-state — devices. Electron 

dynamics. Vacuum devices. 

An introductory book on electronic devices 
it endeavours to present the basic ideas of 
both vacuum and solid-state devices. 
^ Many worked examples from past 
examination papers are provided to illustrate 
the application of the fundamental theory. 


Electromagnetism and its Applications—An 

Introduction 

B. Botton (University of Bath). Van 
Nostrand Reinhold, Wokingham, 1980. 
21 x 29 cm. 156 pages. £10 (clothbound) 
£4.50 (paperback). 

CONTENTS: Static electricity or when to keep 

your shirt on. Charges and forces. Electro- 

static copiers and chimneys that don't smoke. 

Electric field strength and potential energy. 

Electric potential difference. Lightning. kugel- 

blitz and when to stand on one leg. Gauss’ 


law. Motors and rockets. The effect of 
dielectric materials on the electrostatic field. 
Tracking, treeing and bushings. Capacitance. 
Charge-coupled devices. Energy in electro- 
static systems. Electrostatic generators. Diag- 
nostic test. The magnetic field. Navigators all 
at sea. Magnetic flux density. The free-space 
constants. The effect of iron on the steady 
magnetic field. Bubble and squeak. Ampere's 
circuital law. This is your local displacement 
current —here is the news. Magnetic circuits. 
Some small machines. Permanent magnets. 
Do. do, do, do, do you remember? The Biot- 
Savart law. Electromagnetic induction. The 
new ice age. Mutual and self inductance. 
Energy storage in the magnetic field. 
Machinery. The transformer. More about the 
transformer. Hysteresis and eddy-current loss. 
And finally, Maxwell. 

This textbook covering the basic principles 


i of electromagnetism is for first year under- 


graduates and students on diploma courses. 
The chapter headings. quoted above, indicate 
the informal approach. 


Radar Anti-jamming Techniques 
Edited by M. V. Maksimov. Artech House 
Dedham, MA 1979. 15.5 x 23.5 cm. 421 

pages. £29. 
CONTENTS: Natural radio interference Man- 
made radio interference Mutual interference 
and electromagnetic compatibility of 
electronic systems. General description of 
radio interference control techniques Pro- 
tection of receivers from overloads and 
cancellation of radio interference . Spatial, 
polarization, frequency and phase selection. 
Time and amplitude selection. Functional, 
structural and combined selection. Combined 
coordinate sensors. 

Translated from the Russian, this book by 
a group of Soviet electronic engineers is 
intended primarily for those engaged in the 
development and operation of various kinds 
of electronic equipment. It may also serve as a 
textbook for students of advanced radio 
engineering courses. Methods of protecting 
radio systems from radio interference are 
described in detail and wide attention is 
devoted to methods of protecting systems 
from jamming. 


Handbook of Microcircuit Design and Appli- 

cation 

DaviD F. Stout (Ford Aerospace & Communi- 
cations Corporation). McGraw-Hill, New 


York, 1980. 15.5x24 cm. 499 pages. 
£18.20. 
CONTENTS: Introduction to  microcircuits. 


Digital families. Combinational circuits. Flip- 
flops and shift registers. Read-write random- 
access memories, Read-only memory devices. 
Data-word sorting and checking. Sequential 
circuits. Digital signal generation. Serial 
communications — microcircuits. Parallel 
communications microcircuits. BCD circuits. 
The microprocessor and the microcomputer. 
Microprocessor controlled A/D converters. 
Microcomputer-based traffic systems. Key- 
board scanner circuits. Application for the 
6800 microprocessor: A polyphonic music 
synthesizer. Special purpose LSI circuits. 
Operational amplifiers. Active filters. 
Regulators. Nonlinear analog microcircuits. 
Analog signal generation. Sampling and 
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multiplexing ^ circuits. | Analog-to-digital 
converters. — Digital-to-analog converters. 
Phase-locked-loop circuits. 

A comprehensive handbook to help solve 
problems relating to the design and 
application of most classes of digital and 
analogue circuits. 


Radar Electronic Counter-countermeasures 
STEPHEN L. JOHNSTON (Editor). Artech 
House, Dedham, MA, 1980. 21 x 27.5 cm, 
546 pages. £29. 
CONTENTS: Introduction to radar electronic 
counter-countermeasure. General ECCM. 
ECCM in the transmitter. ECCM in the 
antenna. ECCM in the Antenna: Sidelobe 
blanker/canceller. ECCM in the receiver 
(general. ECCM in the receiver signal 
processor. ECCM in the system, operational, 
and other ECCMs. ECM/CCM, ECCM 
efficacy, simulation and analysis. 

These 52 reprints cover a time span from 
1946 to 1978. About one-third are by non-US 
authors from Chile, Denmark, Germany, 
India, Italy, Japan, Netherlands, Sweden and 
including several from the United Kingdom. 
Five papers were originally presented at the 
1978 Military Microwaves conference in 
London. 


Electrical interference in Electronic Systems: 
Its Avoidance within High-voltage Substations 
and Elsewhere 
R. E. MARTIN (Consultant). Research Studies 
Press, Forest Grove, OR, 1979. 14.5 x 22 
em. 198 pages. £11.30. 
CowrENTS: Interference ingress. Interference 
measurements in substations. Susceptibility of 
computers to interference. Susceptibility of 
solid-state circuits to noise and to damage. 
Typical values of attenuation required. Check 
list of desirable precautions. Test speci- 
fications and methods. 

Intended to be used in conjunction with 
other documents to facilitate interference 
immunization by isolation, screening, filter- 
ing, earthing and layout design for electronic 
systems in electrically noisy environments. It 
deals particularly with the use of electronic 
equipment in high-voltage electricity sub- 
stations. and outlines the way in which 
harmful interference may enter the electronic 
systems as well as indicating appropriate 
measures for the control of the effect of such 
interference. It is not concerned with inter- 
ference effects external to the station. 


A Retrospective Technology Assessment: 

Submarine Telegraphy—The Transatlantic 

Cable of 1866 

V. T. Coates and B. FINN (National Museum 
of History and Technology). San Francisco 
Press, (547 Howard Street, San Francisco 
94105 USA.) 1979. 17 x 25 cm. 264 pages. 
$8.50. 

A methodological study on a remarkable 
technological development of its time namely 
the introduction of submarine cables with 
particular reference to the transatlantic cables 
of the 1860s. Submarine cables are shown to 
have brought about significant changes in 
many institutions ranging from shipping and 
commerce to diplomacy and military organ- 
izations. 
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The Engineering of Microprocessor Systems. 

Guidelines on System Development 

Electrical Research Association, Leatherhead, 
Surrey. Pergamon, Oxford, 1979, 21 x 29.5 
cm. 175 pages. £7.00. 

ConTENTS: A guide to the jargon. The 

significance of microprocessors for industry. 

Choosing the right technology. Selecting the 

microprocessor. The development process. 

Setting up the development laboratory. 

The book is primarily for managers and 
other users who have the responsibility for 
microprocessor system developments but who 
may lack direct experience in this field. 


Handbook for Radio Engineering Managers 
J. F. Ross (Registered Professional Engineer 
in Australia). Butterworths, London, 1980. 
16 x 24 cm. 947 pages. £35. 
ConTENTS: Management and organization. 
Engineering economy. Safety practice. Fires 
in radio installations. Environmental aspects 
in radio engineering. Specification and 
contract administration, 

Designed to help anyone associated with 
the management of radio projects, services 
and facilities from the large complex radio 
stations and networks down to the small 
station. It describes managerial factors and 
aspects associated with the design, 
installation, commissioning, operation and 
maintenance of radio engineering services and 
facilities with particular emphasis on 
efficiency, productivity, budgeting, 
organization, safety practice and environ- 
mental obligations. 


The VDT Manual 
A. CaKIR (Technical University of Berlin), 
D. J. HART (/nca-Fiej Research 
Association) and T. F. M. Stewart (Uni- 
versity of Loughborough). IFRA, 
Darmstadt. 1979. 2] x 29.5 cm. 250 pages. 
ConTENTs: VDT basics. Light, vision and the 
Optical characteristics of visual displays. 
Ergonomic requirements for vdts. Ergonomic 
requirements for vdt workplaces. The health. 
safety and organizational aspects of working 
with vdts. 

Despite the vast amount of literature on 
computers, little has been written about the 
people who use them and the conditions 
under which they are used. It sets out to 
analyse the ergonomic, health and safety 
aspects of working with visual display 
terminals. 


Experimenter’s Guide to Solid State 
Electronics Projects 
ALFRED W. Barper. Parker Publishing 


Company, New York, 1980. 16x23 cm. 

214 pages. £9.70. 
ConTENTS: Experimenting with crystal 
detectors, diodes and rectifiers. Projects and 
experiments with transistors. Building and 
experimenting with solid state power supplies. 
Experimenting with silicon controlled recti- 
fiers and triacs. Building photo detectors, 
photo couplers and optical fibre devices. 
Experimenting with solar cells. Experimenting 
with field effect transistors. Project: build 
your own experimental laboratory, Experi- 
menting with linear integrated circuits. 

A collection of easy-to-follow electronics 
projects and tests for experiments with solid- 


State devices such as solar cells, fibre optics, 
l.e.d.s and both linear and digital i.c.s. 


The Art of Electronics 
PauL Horowitz (Harvard University) and 
WINFIELD HILL (Sea Data Corporation). 
Cambridge University Press, 1980. 
17.5x25 cm. 716 pages. £35 (Hard 
covers), £12.50 (Paperback) 
ConTENTS: Transistors. Feedback and oper- 
ational amplifiers. Active filters and osci- 
llators. Voltage regulators and power circuits. 
Field-effect transistors. Precision circuits and 
low-noise techniques. Digital electronics. 
Digital meets analog. Minicomputers. Micro- 
processors. Electronic construction  tech- 
niques. High-frequency and high-speed tech- 
niques. Measurements and signal processing. 
A text/reference book that emphasizes 
electronic circuit design techniques and 
scientific measurements. It begins at a level 
suitable for those with no previous experience 
in electronics and takes the reader through to 
a reasonable level of design proficiency. 


The Microelectronics Revolution 

Tom Forester (Editor). Basil Blackwell, 
Oxford, 1980. 15.5 x 23.5 cm. 589 pages. 

ConTENTs: The microelectronics revolution. 

Economic and social implications. The micro- 

electronic age. 

The aim of the book is to bring together 
natural and social scientists in a positive way 
on a uniquely important topic. The contri- 
butors, all experts in their fields, explain in 
turn the origins and nature of micro- 
electronics and the increasing use of micro- 
processors in everyday products. It details the 
impact of the new technology on society, the 
consequences for employment and the impli- 
cations for industrial relations. 
Broadcasting—An Introduction 
Jonn R. BITTNER (DePauw University). Rustic 

Hall, Englewood Cliffs NJ, 1980. 18 x 23.5 

cm. 508 pages. £10.35 (cloth). 
CONTENTS: Broadcasting as mass communi- 
cation. The beginning of wireless. The 
development of modern radio and television. 
Radio waves and the spectrum. Broad- 
casting’s use of microwaves and satellites. 
Cable. Networks and syndication. 
Educational and public broadcasting. Tele- 
vision in business and industry. International 
broadcasting. Early attempts at government 
control. The Federal Communications Com- 
mission. Control of broadcast programming. 
Regulation station operations. Cable, 
satellites and future perspectives. Advertising 
and economics. Broadcast ratings. The 
research process. The broadcast audience: 
approaches to studying effects. Using the 
library to learn about broadcasting. 

The book is designed for use in single term 
introductory broadcasting courses in the US, 
being directed toward students who aspire to 
be either ‘responsible consumers of broad- 
casting in society’ or practising professionals 
in broadcasting or related fields. The book 
covers mainly the Northern American 
continent but it will provide interesting 
reading for the British student. This is despite 
the fact that the only reference to J. Logie 
Baird is in connection with his experiments in 
video recording and it completely ignores the 
BBC's broadcasting of television in 1936, and 
indeed until after the War! 
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Microprocessor-based 
digital filters 
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SUMMARY 


The concepts of microprocessor-based digital filter 
implementation are presented via a case study of a low- 


pass filter designed using the bilinear z-transform method. 


Implementation aspects associated with an 8-bit general- 
purpose microprocessor-based realization and a single- 
chip signal processor are discussed and resuits included. 


* Systems and Instrumentation Division, Preston Polytechnic, 
Preston, PR1 2TQ 
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1 Introduction 
A digital filter may be thought of as being simply the 
equivalent of an analogue filter implemented using a 
digital processor. It may, however, have a number of 
advantages in many applications, namely: 

It does not drift. 

It can handle low frequency signals. 

Frequency response characteristics can be made to 

approximate closely to the ideal. 

It can be made to have no insertion loss. 

Linear phase characteristics are possible. 

Adaptive filtering may be achieved. 

Digital word-length may be specified by the designer. 


The digital processor used for filter implementation 
may take the form of hardwired dedicated hardware, or a 
minicomputer system, or a  microprocessor-based 
system. The latter form has relative merits of low cost, 
small physical size and programmability, and 
consequently discussion of practical implementation will 
be restricted to microprocessor-based systems. 

Realization of a digital filter requires an algorithm 
(linear difference equation) for calculating the discrete- 
time values used for the filtering process. There are a 
number of established design methods for deriving the 
required linear difference equation via the corresponding 
z-plane representation of the digital filter transfer 
function.! ^ The approach may involve the selection or 
derivation of a suitable transfer function, G(s), 
corresponding to an analogue filter which meets the 
appropriate specification. The digital filter transfer 
function G(z) is then derived from knowledge of G(s), 
using an s-plane to --plane mapping. 

In some digital filter designs an objective may be to 
achieve a time-domain specification, in which case the 
time response of the analogue filter is used as the basis 
for deriving G(z). This method is suitable if the digital 
filter is required to have the same impulse-response or 
step-response as the prototype analogue filter G(s). That 
is, the s-plane to z-plane mapping is achieved via the 
corresponding impulse invariant or step invariant? 
design methods. A point to note is that for the step 
invariant design method. if the analogue prototype filter 
has fast rise-time and short settling time, these 
characteristics will be preserved in the digital filter. It 
should also be noted that for both methods, the 
frequency response characteristic of G(z) may not match 
that of G(s), the closeness of the match depends upon the 
frequency response of G(s) being adequately 
bandlimited. 

However filters are often specified in terms of 
frequency response characteristics and in many digital 
filter designs these are the main considerations. 
Consequently time-invariant design methods are 
generally not satisfactory in such cases. 

An s-plane to 2-plane mapping. which is 
characteristically bandlimiting in its action, is the 
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bilinear z-transform.!? This method is easy to apply 
because it is simply a matter of substituting a function of 
z for each Laplace operator s appearing in G(s). This is a 
common approach to digital filter design and G(s) may 
be based on prototype Butterworth, Chebyschev, elliptic 
or Bessel analogue filters which have the desired 
frequency-domain characteristics. 

It is therefore appropriate to introduce the concepts of 
digital filter design via an illustrative example involving a 
2nd-order Butterworth filter, which is mapped to the z- 
plane using the bilinear z-transform method. 


2 illustrative Design Example 
The transfer function of the normalized 2nd-order 
Butterworth low-pass filter is given by equation (1): 


G(s) = T - 
s^T2s-l 


(1) 

Let the specification for the digital filter be: low-pass 
digital filter cut-off frequency, f, = 100 Hz, sampling 
period, T = 1-6 ms. (f, = 1/T = 625 Hz). The transfer 
function given by equation (1) may be denormalized 
using the relationships given in Table 1. In this example 
we replace s by s/jw,, in equation (1), where œ is the 
radian cut-off frequency of G(s), thus yielding 


| 
G(s) = 
ce) 
+ 
Dea 


The transfer function given by equation (2) may now 
be mapped from the s-plane to the z-plane using the 
bilinear z-transformation. However a non-linear 
frequency distortion is produced by the bilinear z- 
transform. This warping must be taken into account in 
the design process by pre-warping the s-plane frequency 
scale. To achieve this the cut-off frequency of the 
denormalized prototype analogue filter is calculated 
using equation (3) 


2 T 
Og = T tan (2s J 


where «c, is the radian pre-warped cut-off frequency of 


= (2) 
VENE 


Wea 


(3) 


| 472243-84 
z-| 

fis) 
PES 


which simplifies to 


G(z) = 
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-; : 
| + fons x 1250( 


G(z) = 


Table 1 
Denormalizing transformations 


To transform from Substitute for s 


normalized low-pass to: 


low-pass S/a 

EE OgalS 
bandstop Seay ~ Oca VG? + ear aul 
bandpass (5? ea Ocau Sus OL) 


Oy = upper transition radian frequency, Wa = lower transition 
radian frequency 


the denormalized analogue filter, œg is the specified 
radian cut-off frequency of the digital filter, and T is the 
sampling period in seconds. Substituting the specified 
parameters in equation (3) yields 


16x 107-3 


+ tan (200 x ^ ) = 687-2 rad/s 


2 
W 
107? 


- 


ca — 1 6x 
The next step in the design is to substitute this value of 


c, in equation (2), thereby producing the pre-warped 
transformed transfer function, G(s),,,. thus 


G(s) 


WE /2s 

yz 
E e] 
(az) “epee 


472243:84 


Gish = 524 071-85 5447224384 5 
For the bilinear z-transform we use the substitution 
- 2 (54) 
tf Wega 
and for the specified value of T this gives 
zl 
s — 1250 (Et) (5) 


Substituting equation (5) in equation (4) gives G(z), i.e. 


s jazaa 
zt+l (6) 
zix2z41 


68871 — 4-62: + 1-74 
Having derived the transfer function, G(z), for the 
digital filter, it is necessary to establish the location of its 


poles (roots of the denominator polynomial: 
equation (6)), and for stability, ensure that they lie within 
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the unit circle in the z-plane. In this example the poles of 
G(z) are a complex conjugate pair at z = 0:336+j0-374, 
these lie within the unit circle, thereby yielding a stable 
filter. 

It is also appropriate at this stage to check the filter 
frequency response characteristic, G(e""). This is 
achieved by substituting e/"" for z in equation (6) giving 


G(e"") = 


The magnitude |G(e"") and phase /G(e"") 
characteristics were evaluated at several frequencies in 
the range 0 to 400 Hz, and the results are shown in 
Fig. 1. For comparison the frequency response 
characteristics |G(jo))) and /G(j«) of the denormalized 
prototype analogue filter are included. These were 
achieved by substituting j« for s in equation (2). thus 
yielding for a value of f equal to 100 Hz 


3947842 


G(jo) = : 
VO) = (3947842 — o) 4-j8886 o 


The results in Fig. 1 show that up to the filter cut-off 
frequency the digital and analogue filter characteristics 
are closely matched. However for values of frequency 
approaching f,/2 Hz there are significant errors, and 
indeed beyond f./2 aliasing occurs. The characteristics of 
G(z)are dependent upon the sampling period, T, and the 


(a) 


MAGNITUDE (dB) 
8 


aS eee Pon 
100 150 200 250 300 350 400 Hz 


0 50 


ib) 


PHASE (DEGREES) 
S 
o 


G(s) 
zi 
SH G (2) 
-200 
Z2 
Oo 50 100 150 200 250 300 350 400 Hz 


Fig. 1. (a) Amplitude/frequency response for G(s) and G(z). 
(b) Phase/frequency response for G(s) and G(z). 
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(cos 2-0 T +2 cos mT - 1) - j(sin 2:0 T +2 sin wT) 
(688 cos 2wT — 4:62 cos wT + 1-74) - (6:88 sin 20 T —4-62 sin wT) 


matching of G(z) to G(s) may be improved by decreasing 
T if practicable. The direct form of realization of an nth- 
order digital filter may be expressed as equation (8): 


n A YE 
G(z) = je = x 
1+ b! 7 
i1 
where m > l for realizability (8) 


Since z ^ implies a time delay equal to d sampling 


instants, it follows that the corresponding linear 
difference equation may be written as 
I 


y()T = Y aix(n - i)T — Y, by(n—i)T (9) 
i-0 i=l 
Converting equation (6) to the corresponding linear 


difference equation form, gives equation (10): 


y(1)T = 0145x(n)T +0291 x(n — 1)T -0-145x(n — 2)T + 
*0671y(n— 1)7 —0253y(n — 2)T (10) 


In Section 3 two forms of microprocessor system 
implementation of the filter are discussed. 


3 Microprocessor Implementation 

The processor word length required to maintain 
stability, assuming truncation rather than rounding of 
the filter coefficients, may be determined using the 
method reported by Kuo and Kaiser.! ! 

They showed that for an nth-order filter having 
distinct poles at (cos «,T —j sin «y, T), where k is an 
integer in the range n > k > 1, the number of processor 
bits, w, must be 


w — smallest integer exceeding 


See 
f -iog EZ II en 
k-1 


Evaluating equation (11) for the complex conjugate 
poles obtained from equation (6), yields w = 2. This is 
not a surprising result considering that the poles are 
located well away from the circumference (stability 
boundary) of the unit circle. There is no guarantee that 
with this value of w the filter will meet its specified 
frequency response characteristic. Generally, a rule-of- 
thumb is employed whereby normally four or five bits 
are added to w in order to take account of quantization 
and computional errors. This suggests that for this 
example it may be appropriate to consider 
implementation using an 8-bit wordlength. 


(11) 
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MPU 


Input Signal 
x(t) 


T 


Start Conversion 


2 Phase 
Clock 


(M6800) 
PIA 
A/D jux 
^ 
CA 


EPROM 
(2708) 


Output Signal 


2 


Convert Command Signal 


Status Signal 
0 


Conversion Complete 


Fig. 2. Block diagram of 8-bit microprocessor system. 


A typical 8-bit microprocessor-based digital filter 
configuration is shown schematically in Fig.2. The 
flowchart for the filter implementation is shown in Fig. 3. 

In computing r(n) T (eqn. (9)) the product terms in the 
equation involve evaluating the multiplication of two 8- 
bit binary numbers, thereby producing a 16-bit result 
which would normally be rounded or truncated and 
represented by the most significant 8-bits. This obviously 
introduces computational errors which may degrade the 
filter performance. In the design example, discussed in 
Section 2, with 8-bit representation, the rounding or 
truncation effect is insignificant because the poles of the 
filter are located well away from the circumference of the 
unit circle, and a considerable change in the 
representation of the product terms would be required to 
cause a significant degradation of the filter frequency 
characteristics. 

The multiplication may be achieved using a software 
subroutine or a hardware multiplier, the latter being 
generally faster in operation when compared with the 
software subroutine. However. today the cost/bit of 
r.o.m. (or e.p.r.o.m.) is not a serious limitation to the 
digital filter designer, and consequently the look-up table 
method of implementation has the advantage of 
relatively fast operating speed, and it is easily 
programmed. 

The practical results shown in Figs.4 and 5 were 
obtained from an M6800-based realization of the 
second-order low-pass digital filter. Figure 4 shows the 
digital filter waveforms and Fig. 5 is a spectral analysis of 


426 


the output signal, which illustrates the 
characteristic and aliasing effects. 

In terms of chip count the microprocessor-based 
digital filter system shown in Fig. 2 may be reduced by 
replacing the 6800 m.p.u., 6810 r.a.m. and 2-phase clock 
circuit by a 6802 m.p.u., which has an integral clock and 
r.a.m. In contrast, for some practical applications, the 
system may need expanding to provide a multiplexed 
input and/or a multiplexed output. 

An alternative form of realization may be achieved 
using the single-chip Intel 2920 signal processor.!? This 
n.m.o.s. chip has ‘on board’ sample-and-hold circuits and 
d/a and a/d converters. It has its own e.p.r.o.m., scratch 
pad memory and special instruction set for signal 
processing. The functional block diagram of the 2920 is 
shown in Fig. 6. 

Under program control, one of four possible inputs is 
selected and then sampled and held, and the signal is 
then converted linearly to a digital word with up to 9-bits 
of resolution (includes sign-bit). An internal digital to 
analogue register (DAR) accumulates the digital word 
until the conversion is complete, giving the required 
x(n)T. The word is then loaded into the scratch pad 
memory to be used in the subsequent evaluation of 
yG1)T. The DAR is also used to drive the d/a converter to 
any of eight analogue outputs via the output 
demultiplexer and sample-and-hold amplifiers. An 
interesting feature is that the analogue and digital 
operations can be executed simultaneously, and 
consequently this results in a relatively short instruction 


low-pass 
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Clear: 
x(n-1)T, x(n-2)T, 
y(n-1)T, y(n-2)T 


Convert 
Signal 


Generate 
Command 


Wait 
for interrupt 
(A/D Complete?) 


Fig. 3. Flowchart for M6800-based digital filter. 


set, particularly suited to digital signal applications. 

It should be noted that in comparison with the 8-bit 
microprocessor-based filter having a sampling time of 
1-6 ms, the Intel 2920 can sample approximately ten 
times faster, thereby increasing the baseband of the input 
signal. 

A 2920 realization of the low-pass digital filter, with 
T =0-192ms, has been achieved and the results 
compare favourably with those shown in Figs. 4 and 5. It 
is outside the scope of this paper to present details of the 
program development and the corresponding assembly 
language listing—however readers requiring further 
information are encouraged to contact the authors. 
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Fig. 4. Waveforms for M6800-based filter: 
(a) input signal, x(r): 
(b) output signal from d/a converter, y*(r); and 
(c) output from reconstruction filter, y(t). 


Fig. 5. Results of spectrum analysis carried out on the 2nd-order 
Butterworth low-pass filter 
Spectrum analyser settings: 
vertical scale— 10 dB,cm 
horizontal scale— 100 Hz/cm 
bandwidth — 10 Hz. 


4 Concluding Remarks 

This paper has attempted, by using a relatively simple 
case study, to illustrate the main characteristics of a 
digital filter implementation using microprocessor-based 
systems. It has been shown that a general-purpose 8-bit 


processor, with the necessary peripheral signal 
conditioning elements. operated satisfactorily. In 
comparison the single-chip 2920 signal processor 


implements the desired filtering action, but has an added 
advantage of requiring less hardware. Furthermore the 
2920 can easily deal with multiplexed input signals, and 
offer demultiplexed output signals. Undoubtedly the 
2920 offers designers an opportunity to develop single- 
chip signal processing systems. 
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SUMMARY 


This paper defines a region of operation, calling it the 
pseudo-linear region, where a phase-locked loop (p.1.l.), 
behaves essentially in a linear manner. The results obtained 
for a loop with an active second-order filter suggest a 
very simple empirical definition for this region as a 
function of the parameters of the loop (č and «,), and 
classifies the regions of operation of the p.l.l., linear, 
pseudo-linear, non-linear and cycle slipping. 
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1 introduction 

The phase-locked loop (p.1.1.) which is incorporated in 
one way or another in most up-to-date communication 
systems has been extensively studied, either in the region 
where it is in the process of acquiring lock and is 
therefore cycle slipping,' or as a linear system.” 

In actual fact the linear analysis is not really valid 
except for very small phase error, say less than 5 degrees, 
putting a severe limitation on the usefulness of such 
analysis. This leaves out what probably is the most 
important region of operation of the p.l.]., i.e. when it is 
in lock, but having large phase errors. This paper defines 
a region, calling it the pseudo-linear region, where the 
behaviour of the loop, although non-linear, approxi- 
mates closely to the linear model. The definition is based 
on one of the parameters of the transient response of the 
loop, namely the settling time, which in most 
applications is the parameter of greatest interest. If the 
frequency of the input to the loop is suddenly changed, 
the frequency of the output will try to track the input. 
The settling time is that necessary for the output 
frequency to reach, and remain within, +2% of the 
frequency change from the frequency of the input. For 
small changes in frequency, the loop is linear, and the 
settling time is constant. For greater steps in frequency, 
the loop exhibits non-linear behaviour, and the settling 
time increases until a certain step value reached or 
exceeded, when the loop loses lock and begins to cycle 
slip and the settling time increases very sharply. (This is 
shown clearly in Fig. 6.) The pseudo-linear region is 
defined as that region where the settling time does not 
increase by more than 10% of its value for the linear 
model. This puts an upper bound on the step of 
frequency which may be applied (360 rad/s in Fig. 6). The 
pseudo-linear region was studied for a second-order loop 
with active filter which is the optimum estimate? for a 
step of frequency, and is thus most widely used, together 
with the second-order loop with the passive filter, which 
anyway approximates quite closely to it.* Due to the 
non-linear nature of the problem, an exact analysis is not 
possible. The response of the non-linear loop was solved 
on a digital computer, using the Runge-Kutta method, 
and the pseudo-linear region is described in terms of the 
natural frequency («,) and the damping ratio (č). A very 
simple empirical relationship was obtained relating the 
frequency step with œw, and £. This enables the designer 
to choose his parameters so as to stay in the pseudo- 
linear region, where the loop may be considered 
approximately linear. These results tie in with an 
empirical formula recently given by Gardner? for the 
frequency step at which cycle slipping begins, and thus 
the regions of operation of the p.l.l. are clearly classified 
as shown in Fig. 9. 

2 Theoretical Background 
The block diagram of a p.1.l. in lock is shown in Fig. 1, 


where: 


Input filt) = A, sin (ot +9;) 
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Fig. Il. Block diagram of the phase-locked loop circuit. 


Output 


Jt) = Az sin (or + ; + e.) 


and 
9 = Pi- Po. 


The filter shown is active and the loop is second 
order*, i.e. the d.c. gain of the loop is theoretically 
infinite. 

The equation describing the system is then 


@+(K cos q)ó + Ka, sin ọ = ĝi. (1) 


This is a non-linear differential equation in o. If ọ is 
very small, then it can be reduced to 


$+Ko+ Kaze = $i. (2) 


2.1 Linear Analysis 

Equation (2) describes the linear model of the p.l.l., 
where the output of the phase detector in Fig. 1 becomes 
K,q@ instead of K, sing and its solution in the time 
domain for a unit step input in frequency can be shown 
to be: 


e Ont 


Aw = rus (ó cos ót — £c, sin dt) (3) 


where Aw is the error frequency between input and 


output. 
wna = ./Ka, 


ô= wn /1- E. 


The root locus can be very useful in the linear analysis 
of a p.1.1.* The open-loop transfer function of the loop is 
given by K(s+a,)/s? and the root locus is shown in 
Fig. 2, which is a circle of centre —a,, and radius a,. The 
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Fig. 2. Root locus plot of the p.l.l. with active second-order filter. 


value of K determines the position of the closed-loop 
poles: 


K — = 
St iy aKa Ke (4) 
--—atjb, 
where 
K e 
Ho b = 3,/4a,K —K? 
and 
o, = a? +b? = J/Ka; 
K 
é => sin B = Es = — — 
wo, 20, 


The settling time 7; is the time necessary for Aw in 
equation (3) to remain within - 0-02. Since it is quite 
laborious to obtain T; from equation (3), quite frequently 
it is evaluated with a good degree of approximation (for 
é = 0-707) as: 

T=- 
a 

Figure 3 is a plot of T, versus č for w, = 100,/2 rad/s 

for the exact value of J. 


2.2 Non-linear Analysis 
If the frequency step applied is large, the linear analysis 
given above no longer holds, and a solution of equation 


theoretical 
xxx practical 


wa = 100 /2 rad/s 


Fig. 3. Plot of settling time versus damping factor (č) for the linear p.1.1. 
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Fig. 4. Trajectories in the phase plane for the p.l.l. 


(1) becomes necessary. Since it is a non-linear differential 
equation, an exact solution is not possible. Trajectories 
in the phase plane" can be obtained, and Fig. 4 shows 
such a trajectory. However, it is not easy to obtain the 
settling time from such a trajectory, since time is an 
implicit variable: in the phase plane. To obtain a 
definition for the pseudo-linear region, equation (2) was 
solved on an IBM 1130 digital computer, using the 
Runge-Kutta method, for different values of input step 
in frequency, and for different values of € and w, (and 
thus different values of K and a,). From the solution, it 
was possible to obtain 7; as defined in Section 2.1. 
Figure 5 shows a plot of 7; versus input frequency step 
value for one value of w,. It is a family of curves, each 
curve for one value of č. Similar data were obtained for 

wa = 50, 50,/2, 100, 1004/2, 200, 200,/2, 

and 400 rad/s. 
3 Pseudo-linear Region 
The pseudo-linear region denotes that region where the 
behaviour of the p.1.1., although not exactly linear, can be 
approximated to that of the linear model. Since non- 
linearity exhibits itself in the transient response of the 
loop, one of the parameters* of the transient response 
must be chosen for a measure of linearity. The settling 
time is the obvious parameter to choose since it is the 
parameter normally of greatest interest. It is expected 
that similar results will be obtained if any of the other 
parameters are chosen, e.g. peak time or maximum per 
unit overshoot. 

The loop is defined as pseudo-linear when its settling 
time is not greater by more than 10% than the settling 
time of the linear model. This requires that the input 
frequency step should not exceed a certain value y, which 
is then a function of č and c. Figure 6 shows a plot of T; 
versus the step Aw, when w, = 100, /2 rad/s, č = 1/2 
and shows how y is evaluated. 
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3.1 Empirical Relationship 

As shown in Section 2.2, families of curves similar to 
Fig. 5 were obtained for different values of w,, and for 
each curve y was evaluated. This enables Fig. 7 to be 
plotted, giving y versus č for different values of c, from 
which constant y contours are obtained as shown in 
Fig. 8. Here, the pseudo-linear region is well defined 


F-o1 


0:2 


wp = 100/2 rad/s 


Ts (s) 


0:5 
0:707 


09 12 1-414 J1-7 2 


y (rad/s) 


Fig. 5. Plot of the settling time versus the step of input frequency. 
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Fig. 7. Plot of y vs. č. 


and, for a certain value of y, it is the area bounded by the 
constant y curve away from the origin. 

The hyperbolae of Fig. 8 indicate that if y is 
normalized with respect to w, and averaged for all values 
of č, it may be possible to eliminate c, and obtain: 


Pues ( 4 ) scr 
On averaged over ¢ 


When this was carried out, it was found that y/w, was 
fairly constant for any one value of £, and hence it was 
possible to average it to obtain y, as a function of č only. 


(5) 


This was plotted in Fig.9, yielding the linear 
relationship: 

y, 124196 (6) 
Or 


y = (12419£)o,. (7) 


Due to the non-linear differential equation describing 
the phase-locked loop, it is not possible to justify in an 


| 
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L ———M— er 
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Fig. 8. Constant y contours. 
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exact manner the relationship given by equations (6) or 
(7), and as such it is proposed as an empirical 
relationship valid for the region investigated and shown 
by the full line in Fig. 9. However, it can be argued 
intuitively that the relationship is quite reasonable. For 
very small values of the damping ratio č, the step in 
frequency which the loop can accept, while remaining 
linear, is directly proportional to the bandwidth of the 
loop and hence to w,. This is the first term of equation 
(6), which gives the constant of proportionality. As the 
damping ratio increases, the loop will accept a larger step 
in frequency, and will have a larger bandwidth. This 
effectively means that the constant of proportionality 
increases by a factor proportional to the damping ratio 
as given by equation (6). 

It should be noted that since č is defined as K/2o,, it 
can have a value greater than unity, and in fact the 
results given are up to a value of € = 2. 


Fig. 9. Operating regions of the p.l.1. 


Equations (6) or (7) may now be taken as a definition 
of the pseudo-linear region, i.e. for any value of £ and w,, 
y from equation (7) gives the upper bound for the 
step in input frequency which may be applied, for the 
loop to behave in essentially a linear fashion. For 
(c 152. which is a very common value to use, 
equation (7) states that the step in input frequency 
should not exceed about 2:5w,. The trajectory corre- 
sponding to this is shown by a heavy line in Fig. 4. 


3.2 Operating Regions 
It is now possible to define the various operating regions 
of the p.l.l. in a y, vs. č plane as shown in Fig. 9. The x- 
axis, i.e. y, — 0, is the linear operating region. The area 
bounded by the x-axis and equation (6) is the pseudo- 
linear region, where the loop may be approximated to 
the linear model. 

The empirical relationship given by Gardner,? based 
on Viterbi's results," 


Ya = 1'8(+1) (8) 


gives the lower bound for cycle slipping. The area 
between the cycle slipping and the pseudo-linear region, 
is the non-linear region, where the loop is definitely non- 
linear. Since T, rises quite sharply near the cycle slipping 
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Fig. 10. Circuit diagram of the practical p.l.l. 


value of Ac as shown in Fig. 5, the non-linear region 
fortunately is quite narrow. 


4 Experimental Procedure 

A practical phase-locked loop circuit was constructed as 
shown in Fig. 10 to demonstrate the variation of T; with 
Ac. Both the input source and the voltage controlled 
oscillator (v.c.o.) were synthesized signal generators 
operating at 150 kHz. The step of input frequency was 
introduced as a step function to the difference amplifier 
A2 between the active loop filter and the v.c.o. By 
writing the equation of the loop, or by inspection, it can 
be seen that this step of voltage is in effect a step of input 
frequency. An alternative method is to apply a step to the 
external f.m. input of input source, if available. The input 
to the v.c.o. is taken as the output frequency and 
displayed on a storage oscilloscope, from which 7; may 
be calculated. Figure 11 shows a trace of the output from 
the linear p.l.l. for m, = 50, 100, 200 and 400 rad/s. The 


Fig. 11. Trace showing the transient response of the p.l.l. 
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measured values of T, at c, = 100,/2 rad/s are shown in 
Fig. 3 plotted against K for the linear p.l.., and at 
On = 100,/2 rad/s and č = 1572 are shown in Fig. 6 
plotted against different steps of input frequency. In both 
plets close agreement with theoretical values is evident. 


5 Conclusions 

The pseudo-linear region of operation of the phase- 
locked loop was taken as that where the p.1.1. behaves in 
essentially a similar manner as the linear model. It was 
specified from the transient response of the loop as that 
region of operation where the settling time, after a step in 
input frequency, does not exceed the settling time of the 
linear model by more than 10%. By solving the non- 
linear equation of the p.l.]. on a digital computer, it was 
possible to map this region in an œ, versus č plane for 
different values of input steps, y. It was shown that if 
these input steps are normalized with respect to œn, then 
the normalized value is almost the same for any one 
value of € and, thus it was possible to average the 
normalized steps over č calling it y,, and obtaining an 
empirical expression of y, as a function of £ only, which 
turns out to be a linear relationship. This is used as a 
definition of the upper bound for the pseudo-linear 
region. It becomes, therefore, possible to specify the 
maximum step in input frequency which any p.l.l. can 
accept before the linear analysis breaks down, and how 
the parameters of the loop (č and c) may be modified to 
increase this region. À practical model was constructed, 
and showed agreement with the theoretical and digital 
computer results for the measured value of the settling 
time. The work is limited to the second-order loop with 
active filters, which are the types most extensively used, 
and to which many other types may be approximated. It 
is hoped that it will stimulate further work to specify the 
pseudo-linear region for other types of loops. 
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SUMMARY 


The paper develops a methodology for the automatic test 
program generation (ATPG) for a large analogue circuit 
subject to component drifts. An optimization algorithm 
developed selects the best set of intermediate access 
points for a.c. and d.c. tests of the circuit. 

During ATPG, faults are simulated by varying each 
component over a wide range. Voltages measured at all 
available nodes are thus used in the pre-processing stage 
of node selection. An optimization procedure based on the 
discriminating power of measurements and separability 
between fault signatures is used to select the best set of 
frequencies applicable for all nodes in a.c. testing. A pre- 
amplifier and control amplifier of 56 components is 
chosen as Unit Under Test (UUT) to demonstrate the 
ATPG. The fault isolation scheme used is based on the 
nearest neighbour rule. 

The relationship between diagnosability and number of 
test features is shown to follow the customary Pareto type 
curve. However, four carefully chosen nodes and three 
carefully chosen test frequencies are shown to give an 
adequate level of diagnosability. The whole scheme is 
implemented automatically using a minicomputer 
interfaced to the UUT. 
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Nomenclature 


UUT Unit Under Test 

ATPG Automatic Test Program Generation 

dik measurement for kth class at ith feature; 
also, fault dictionary with i= 1,2,... M 
features; k = 1,2,... L fault cases 

ô threshold for ‘common’ attributes 

MT measure of testability 

PMT: partial measure of testability for d.c. 

Ros) nominal response of the UUT, i= 1,2,... 
M frequencies; j = 1,2, ... N nodes 

T response of faulty system, i— 1,2... M; 

j=1,2...N 

Y measurement vector for an unknown fault 

Y, deviation vector at jth node 

X ijk recognition matrix, i — 1,2... M 
frequencies; j = 1,2... N nodes; 
k — 1,2... L fault cases 

Xm Euclidean normalized matrix of X;;, 
normalized over i 

Di discriminatory power of ith frequency at jth 
node 

d(k,, kj) separability measure for fault cases k, and k, 

C, confidence level 

Cy total confidence level 

G percentage confidence level 

dos optimum frequency set 

IN, test frequency set at node N, 

JE frequency set before optimization 

d, nearest neighbour distance for fault case k 

o standard deviation of random noise 


1 Introduction 

The demand for automatic testing and fault isolation in 
analogue circuits and systems is rapidly increasing in a 
wide range of industries. This is because substantial cost 
and time are being expended in manual testing, location, 
and correction of faults caused by defective material and 
manufacturing processes. Presently, the desperate 
shortage of skilled inspectors and repair technicians gives 
added impetus to the use of ATE to increase workstation 
throughput, product quality, and productivity generally. 

Analogue automatic test program generation (ATPG) 
should, in principle, provide a software algorithm which 
receives the circuit or system description as input, and 
then produces a test program which automatically 
isolates a failed part.! This is a particularly important 
area of research now that packaging densities of 
integrated circuits per board are so high that manual test 
program generation is out of the question. 

In fault isolation, there is a tenuous relationship 
between the provision of access points and increased 
diagnosability (i.e. the ability to correctly isolate the true 
fault). Our particular interest is in maintenance testing 
analogue systems comprising large electro-mechanical 
and electro-hydraulic components as well as electronic 
circuits. In some ways testing such systems raises more 
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Fig. i. Example UUT. 


Table 1 
Nominal component values for UUT 


Pre-Amplifier Control Amplifier 
Resistors Capacitors Resistors Capacitors Potentiometers 
RI S56 kQ R21 1000 C1 220nF R30 2700 CIS 1lp/12V PI 1 k log 
R2 $56 kQ R22 100kQ C2 470pF R31 18 kQ Cló  22nF P2 10 k log 
R3 56 kO R23 1000 — C3 I004F/16 V R32 18kQ CI7  22nF P3 100k lin 
R4 S6kQ R24 33kQ C4 27pF R33 33k 
RS S6 kO R25 27kQ CS 10pF/16V R34 56kQ C18 82nF 
R6 1kQ R26 6800 C6 12nF R35 470kQ C19 82nF — ———————— 
R7 33Q C7 12nF R36 12kQ C20 IOuF/12V Semiconductors 
R8 470 KQ C8 82nF R37 47kQ C21!  ISnF op 
R9 100 kQ C9 220uF/4 V R38 3300 C22 Iy/FI2V TRS 3C179C 
RIO 68 kQ CIO 10pF/25 V R39 5600 C23 IOuF/12V TR6 3C109C 
RI] 15 kQ CI] 10pF/25 V R40 3900 (or 
RI2 68 kQ C12 100 pF/25 V equivalents) 
RI3 22 kQ - — 
RI4 œ Semiconductors 
RIS 2700 
RI6 10 kQ TR!  BCI79C 
RI7 33 kQ TR2 BCI79C 
R18 39 kQ TR3 BCI09C 
RI9 27 kQ TR4  BCIO9C 
R20 56 kQ Z1 (or 
equivalents) 
Zener diode 
18 V 400mW 
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severe problems than circuit testing. For example, the 
provision of additional access points is considerably 
discouraged, on the grounds that the provision of any 
extra components will degrade reliability.? Also, there is 
no guarantee that a randomly chosen access point will 
shed further light on the likely cause of failure. Even a 
well-chosen access point may be of little value for 
diagnosis if the noise level is high. Attempts to reduce 
noise effects by integration will mean a longer test time 
with consequential increased wear and reduced test 
station throughput. There is thus a need for the designer 
to be provided with a means for selecting the minimum 
number, and indeed, the best set of additional access 
points. ATPG should be so engineered as to greatly 
assist this process. 


2 Fault Isolation from a Limited Number of 
Access Points 


Over the past few decades a number of approaches for 
fault isolation in large scale analogue systems from a 
limited number of measurements have been proposed. 
These include the use of both time and frequency domain 
and statistical, non-deterministic and deterministic 
methods.?-? Historical reviews of analogue fault location 
are available elsewhere?'!? and need not be repeated 
here. In this paper a fault diagnosis scheme which 
isolates component drift faults via a functional testing 
approach injecting both d.c. and a.c. signals is used to 
illustrate how access point selection may be improved. 

Catastrophic failures are not considered here. There 
are a number of schools of thought as to which fault 
types are most likely to occur, as design, environmental 
and maintenance crew quality all affect the issue. Our 
experience is that the uncertainty caused in analogue 
systems by the combined effects of measurement noise 
and normal production variances makes drift fault 
isolation the more complex problem. In contrast, their 
catastrophic failures are often of such a nature as to 
rarely require sophisticated test methods. 


3 Problem Statement 
The problem we seek to solve is: ‘given an adequate 
mathematical model describing an analogue system, 
including production tolerances and noise levels, what 
minimum set of test failures should the designer select so 
as to provide an acceptable level of fault diagnosability?' 
In the context of this paper, a test feature is a single 
measurement data point. Thus at two extremes a four 
feature set could be four d.c. measurements made at 
different access points in the system, or gain and phase 
measurements at two selected frequencies based on single 
input-output access only. To solve this problem, the 
designer needs to be provided with ATPG which will 
provide a numerical measure of the *worth' of the feature 
set, so that sufficient confidence can be placed in the test 
design before volume production commences. Additional 
factors to be considered in ATPG include the most 
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relevant stimuli for the particular unit under test (UUT), 
ie. static testing, swept frequency testing, multi- 
frequency testing etc., and the test data processing 
required, both for feature estimation and subsequent 
fault location.!! 

In this paper a fault is arbitrarily defined as a change 
in a component value which results in the system's 
performance lying outside test gates. Therefore, each 
component can generate many faults. This definition 
caters for those system users who believe that an estimate 
of the value of a failed component is an aid to quality 
assurance monitoring. 


4 Contribution of the Present Paper 
The authors feel that the paper will be a useful aid to 
designers, since it is based on experience gained by à 
group of engineers on a range of practical systems. The 
results presented on feature selection using only input- 
output dynamic measurements have been widely tested 
by simulation, supplemented where possible by 
hardware experimentation. There appears to be an 
acceptable degree of correlation between the 
diagnosability achieved in large scale simulations, and 
that predicted a priori from data used to select the 
optimum feature set. 

This paper differs fundamentally from some others in 
the literature in that access points are selected to enhance 
diagnosability, as distinct from necessarily simplifying 
computational procedures. It therefore sheds light on the 
engineering ‘worth’ of an access point. 

D.c. circuits have been less widely studied by the 
authors, but the UUT considered as an example herein 
has been implemented in hardware and software form, 
and considerable experience gained thereon. By choosing 
a relatively common circuit type, problems of military 
and commercial confidentiality have been avoided: 
Although research into test techniques is analogous to 
hitting a moving target, it is thought that the results will 
be sufficiently general to withstand at least the immediate 
effects of galloping obsolescence. ; 

Insofar as the ultimate testing and fault diagnosis of 
hardware is concerned, we have concentrated on 
techniques suitable for mini/microcomputer 
implementation. The same constraint has :not been 
applied at this stage to software developed to provide 
simulation data for access point optimization. This 
inefficiency is accepted, and the approach will be refined 
at a later date. 


5 Unit Under Test 

The circuit selected as the UUT to demonstrate the 
access point selection technique developed later in the 
paper is a pre-amplifier and control amplifier circuit of 
56 components. The circuit is shown in Fig. 1 and the 
component values in Table 1, and it can be used for 
many applications depending upon the input selector 
which preceeds the  pre-amplifier. Together with 
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additional switching inside the feedback loop the latter is 
capable of coping with both high and low level 
signals.'?.!? The circuit can therefore be optionally 
matched to any normal input source by selection of fixed 
resistors and/or equalization networks. 


6 Check-out Tolerances 

In fault diagnosis, we must first define the demarcation 
line between ‘healthy’, i.e. systems in a state of 
operational readiness, and ‘sick’, or failed systems. 
Whereas within digital systems the issue is clear cut (0 or 
1), there is a surprisingly wide band of acceptable 
performance to allow for production tolerances in 
analogue systems.'* In contrast to the usual case of 
assuming parameter independence, we have set check- 
out tolerances at each node using the worst case 
method.'? Circuit analysis is used to determine which 
end of the tolerance range of a component will cause the 
worst performance. 

After the nominal node voltages have been calculated, 
the sensitivities of all node voltages with respect to all 
circuit parameters are estimated. Subsequently (+) 
indicates that worst high performance is achieved by 
selecting the component at the upper tolerance limit, 
whereas ( — ) indicates that worst-case high performance 
is achieved by selecting the component at the lower 
tolerance limit part. Each component value is set at its 
worst-case performance limit, and the equations are 
solved with this unique set of parameters to determine 
worst-case performance. 

A unique set of parameters is required for each 
performance characteristic for the *worst-case high' and 
*worst-case low' situations. This method has been used 
to set up nodal voltage limits due to component 
tolerances typically set at 5%. 


7 D.C. Modelling 

D.c. fault classification concerns only the d.c. level 
setting part of the circuit. Hochwald and Bastian! 
report a method of fault location using only a d.c. 
analysis of electronic circuits. However, they have to 
relate all failures to their d.c. effect and they indicate the 
need to additionally use a.c. analysis for complex 
systems. Current measurements have not been used 
because of the associated complexity of terminal 
disconnections, particularly on sub-networks which have 
to be isolated by special electronic partitioning devices 
.such as those described in Reference 17. 

A 24 V external supply is needed to set the levels in the 
circuit in Fig. 1, while the levels in the control amplifier 
are set by 18 V through the diode Z1. Therefore the 
whole circuit is considered as two independent 
subcircuits, d.c. models of which are shown in Fig. 2 in 
which the nodes are numbered sequentially. This 
partitioning reduces the modelling task. Also fault 
isolation on both circuits can be implemented separately 
if required. 
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Note: 


For component values 
refer to Table 1 


(a) 


(b) 


Fig. 2. D.c. models of UUT. 


8 Selecting Access Points Using the D.C. Model 
It is desirable that the number of access points is kept to 
a minimum so as to reduce the amount of software, 
hardware, test time and cost needed. However, for 
adequate diagnosability, it is not only essential that all 
likely faults in the circuit be propagated to one or more 
of the access points, but also their unique fault 
representation must be maximized. Therefore, it is the 
amount of discriminatory information between the faults 
which is of prime importance in selecting the best set of 
access points. 

Suppose we make a measurement of the ith feature, for 
a system with fault k, and then with fault (k -- 1). Then a 
‘distance measure’ showing the separation between these 
two fault cases is |a,,—a,,41)| where a is a general 
symbol or ‘attribute’ for a measurement of a specific 
condition. If, 


(1) 


we cannot distinguish the two cases and additional 
features are required for diagnosis. The threshold ó is 


lai, — diy « 1) S Ó 
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Form the fault dictionary 
Cn e M measurements 
[ong L fault cases 


Determine the No. of 
common attributes (7) for 


YES 
Calculate 
RMIT. 


All Choose the 
Reset the 
3 i ion access point 
starting combinations set with the 


considered 


ccess point 
a p ? 


highest PMT, 


Fig. 3. Flow chart for the selection of access points for d.c. analysis 


preset and can be varied, based on the data and required 
classification level. Also 6 can be related to measurement 
noise and permitted production variation of system 
components, decreasing with signal-to-noise ratio and 
increasing for large permissible parameter variations. 

A general flow chart for the selection of access points is 
given in Fig. 3. The fault dictionary indicated is simply a 
matrix relating the measurements (M in total), to the 
number of fault cases deemed likely (L in total). 


A ‘quality factor’ or ‘performance index’ based on 
testability is incorporated into the access point selection 
algorithm. Here testability is defined as the characteristic 
of a design which allows the status (operable or 
inoperable) of a system or of its subsystems to be 
confidently determined in a timely fashion.'® This 
definition implies satisfactory determination of probable 
failures and measurements of the likelihood of test 
success (the correct determination of failed status). One 
simple measure of testability (MT) is defined as:'? 


(2) 


MT= [= number of fault cases propagated by Are 


total number of possible fault cases 


However, in the context of fault diagnosis it is necessary 
to maximize unique fault representation thereby 
maximizing the discriminatory information between the 
fault signatures using d.c. and a.c. approaches. So, to 
take account of this, the measure of testability in 
equation (2) is modified and termed partial measure of 
testability (PM T) defined?? by 


PMT is. m 


total number of fault) [ total number of common 
signatures generated fault signatures 


total number of fault signatures generated. (3) 
which is bounded as 0 < PMT,, < 1-0. This measure is 
equally applicable to a.c. analysis. The total measure of 
testability (TMT) is arbitrarily chosen to be the average 
of the two PMT separately obtained for d.c. and a.c. 
testing. 


9 A.C. Analysis 


The a.c. fault isolation scheme for the circuit, the a.c. 
model for which is shown in Fig. 4, is applied after the 
d.c. fault analysis phase and directed at those 
components which are not present in the d.c. analysis. 


= NODES 
= SELECTED NODES 


[m] 
A 


v 


Fig. 4. A.c. model of UUT showing selected access points 
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The a.c. analysis uses the steady state response of the 
circuit subjected to sinusoidal excitations carried out at 
different frequencies over the dynamic range of the 
circuit. 

Although the classical approach of measuring the 
response of a system to completely define its transfer 
function is very time consuming for large scale circuits, 
the a.c. approach, using more than one frequency, does 
contribute information which much reduces the number 
of access points needed. Hence, it is important that the 
number of test frequencies should also be a minimum, 
and chosen for their discriminatory power between 
various fault cases. Conversely, the independence of 
some components to frequency effects inhibits their fault 
location from input-output dynamic tests, and therefore 
the need for additional access points is indicated. 


10 Setting-up the Recognition Matrix 
A recently reported computer-assisted optimization 
technique for feature selection?! is used in the present 
study for the selection of the best set of test features. For 
the same purpose, Tinaztepe and Prywes,?? using a 
program which generates ATLAS coding for the UUT, 
apply a fault approach to optimize the number of tests. 
This program is time consuming, and requires 
engineering expertise in its operation. Although fault 
location in a specific system is frequently possible, using 
only input-output frequency domain measurements,?? 
for the present paper as the multi-test point approach is 
required, the combined set of access points and frequencies 
chosen must be optimized. 
Let the nominal response of the UUT be 
2a. 17 1|2,... M frequencies 
Rromll, j); j21,2,... N nodes 


and the response of a faulty UUT be R;; Then the 
deviation vectors for the jth access point will be 


(4) 


Y = (Rj— Rom) (5) 
For a given fault case k, the deviation vectors are 
Xj, = (Rj, — Roo U)). j= l, 2:5 . .N (6) 


The deviation vectors for all known fault cases 
k—1,2...L are grouped to form a matrix termed 
recognition matrix (X;;,) for 


= 1,2,... M frequencies 
1,2,...N nodes 
k — 1,2, ...L fault cases 


(7) 


~ ~. 
ll 


which is considered as the a priori information of the 
UUT for test frequency selection. (X;,) will be an 
(MN)x L matrix. 


11 Discrimination between Faults 

Optimization in this context is the minimization of the 
number of frequencies and nodes without losing 
discriminatory information between fault signatures. A 
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Lth fault case 


(L~1)th fault case 


Coded deviations Euclidean normalized 


Measurements 


Fig. 5. Calculation of discriminatory information 


Discriminatory information given by: 
D, 2 [hà - h5 9... the R i=1,2,...M 


(see eqn. (8) 
Not all combinations are used 


‘distance’ measure for the discriminatory information of 
the fault signatures in the recognition matrix is one of the 
criteria used in the optimization algorithm and is defined 
as 


L-I 4 
»,-| 5 Oe Xt | fori21,2...M (8) 
oo j=,2,...N 


where Xf, is the Euclidean normalized matrix of X;;, 
over i. Normalization is necessary in order to ensure that 
all quantities are non-dimensionalized and properly 
scaled. Equation (8) does not imply an exhaustive search 
over all possible combinations of faults, frequencies, and 
nodes. It is a much restricted definition which 
considerably reduces computation time during ATPG, 
but which has been found adequate for feature selection. 
Figure 5 is a diagrammatic illustration of the calculation 
of the discriminatory information. 


12 Separability and Confidence Level 

The separability measure used is a distance measure for 
calculating information on the separability of fault cases 
k, and k, at each node as follows: 


M E 
d(k,, kj) = [3 anx | 


i=l 


(9) 


There are L(L — 1)/2 such ‘distances’ for L fault cases at 
each node. For the sake of clarity, the calculation of 
separability measure is illustrated in Fig.6. Due to 
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Coded deviations Euclidean normalized 


Measurements 


Fig. 6. Separability measure between two fault cases K, and K 
Separability measure d(K,, K;) given by: 


M B 
d(K,, Kj) = | (xika — | 
=i 


=[rf+ri+... +r] 
(see eqn. (9)) 
Euclidean normalization, the range of the separability 
measure is bounded and is 
0 < d(k, kı) < 2 (10) 


If for the fault cases k, and k, equation (9) becomes zero, 
the chosen test frequency set is unable to discriminate 
between those two cases. To complete the frequency 
selection algorithm of Fig. 7, a confidence level is defined 
which limits to unity those combinations of fault cases 
with a high separability measure. The confidence level, 


C, = 1 if d(ky, k)2 1 
= d(k,, k,) otherwise (11) 


If C4 is the total confidence level for all combinations of 
fault cases at a node, the percentage confidence level for 
L fault cases is, 


C 
QE l 
[s d us 
Suppose Nj — 1,2,...n nodes are considered and 
fw; represents the optimum set of test frequencies at N; 


node, then the optimum set of test frequencies for all 
nodes is defined for the present study as 


Kop = ENRE : T (13) 


ie. the optimum set is defined as the selected test 
frequencies which are common for all nodes considered. 


13 Test Point Selection in the A.C. Approach 

Having selected the optimum set of test frequencies it is 
necessary to minimize the number of nodes required to 
be monitored for the fault diagnosis scheme. An 
algorithm similar to the d.c. approach is used to select 
the best set of test points. However, the a.c. analysis has 
the additional information from the application of 
various test frequencies and therefore the modification to 
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Recognition matrix for all sampled 
frequencies and available nodes 


M frequencies 


N nodes 
L fault cases 


Eqn. (7) 


7- 


Call rec. matrix at a node [xjx] 


Euclidean normalization [xix] — [xi] 


Calculate discriminatory information Eqn. (8) 


{oit (A), A 2127 £u | xt tina] 


Sum = sume (7) + sume (J +1) 


i (sume (7) -= EM I) 


YES ts 
{Ditt (k), k= 12, 


<sum/p 
ü 


YES NO 
Is 


the same 
p=p+05 


frequency set 
selected before 
? 


Form ‘subset rec.matrix for 
the selected frequencies and 
normalize 


Calculate confidence level [ Eqn. (12) 


is 

maximum 

confidence level 

achieved 
? 


YES 


Choose the set with 
maximum confidence 
level at node / 


NO Are 


eel all nodes 
considered 
2 
YES 


Form the optimum set 
of frequencies for 
all nodes together 


Fig. 7. Flow chart for the optimization of test frequencies 
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Fig. 3 is that nodes are added only after determining 
‘common’ attributes based on all optimum test 
frequencies. 


14 Test Point Selection for the Pre/Control 

Amplifier 
An ECAP software package?^ has been used to calculate 
the node voltages. This package enables the d.c. and a.c. 
nodal voltages of a given circuit to be calculated. The 
user can also increment any parameter and sample from 
a given statistical defined population. Jagodnick and 
Wolfson? have reported a development of ISPICE 
which gives the operator the opportunity of building in a 
constraint on the changes in a particular parameter. 
However, the development is restricted to simulating 
fault nodal voltages and not their selection. Although the 
pre-amplifier and control amplifier are not separable in 
the a.c. analysis, due to the independency of d.c. supply 
voltage requirements they are considered as two 
independent subcircuits for d.c. analysis as shown in 
Figs. 2(a) and (b). In the d.c. approach, the computer 
simulation of both circuits is carried out in parallel. The 
selection of the best set of test nodes is based on the 
sensitivity analysis (fault simulation) of components in 
the circuit. For the present study the fault dictionary so 
formed is a grouping of a coded version of the node 
voltages for known component faults imposed. To 
provide an accurate software/hardware comparison, 
each component value has been varied from —90% to 
+90% of its nominal value in increments of 10%. 

This fault dictionary, although large, would enable 
the designer to consider all likely drift faults. In this 
simulation, each parameter has 18 fault cases resulting in 
18 features, but those drift faults which do not produce 
response deviations and/or whose deviations in response 
do not vary from others are eliminated, which reduces 
some of the redundancy in the matrix. The maximum 
and minimum voltages at each node are noted from the 
matrix and the interval is subdivided and coded 
sequentially. 

It is not essential that the sub-intervals for coding be 
equal but the number of sub-intervals is important. 
Although increasing the number of sub-intervals 
increases the cost of software, selection of small intervals 
provides a better chance of creating unique fault 
signatures. In this study, if the deviation in response is 
near the nominal response, a sub-interval is created and 
coded for every deviation range of approximately 27, of 
the nominal response value. For a larger response 
deviation, a sub-interval is formed only for a deviation 
range of 10% from the nominal response value. This 
variable range assists detection of fault cases causing 
similar deviations to have unique representation due to 
small drifts in component values. 

A coded fault dictionary for the pre-amplifier and 
control amplifier using voltage measurements at 9 and 6 
nodes respectively is formed. The fault dictionary has 79 
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coded fault signatures for the pre-amplifier and 72 for the 
control amplifier. Note that for d.c. analysis the pre- 
amplifier has 8 components considered while the control 
amplifier has 6 components which originally. gave 144 
(8x18) and 108 (6x 18) fault signatures simulated as 
explained earlier. Redundant signatures for each 
component are eliminated to give the reduced number of 
signatures. 

The algorithm described in Fig.3 is 
fault dictionaries and selects various 


applied to the 
sets of nodes. 


(2,3,45,6.7) 
(2,3,4,5,6.7,8) 
(2,3,4,5,6,7,8,9) 


(1,2,45,6,7,8,9) 


(a) 


Nodes 
100 


Percent of unique fault signatures 


80 x 
` — 
70 z) « 
25282 
60 žes 
Sj § 

50 ge~ (b) 


Total no. of fault signatures 72 


| (1,6) 


3 [3 5 6 7 8 9 10 
Nodes 


Fig. 8. Percentage of unique fault signatures achieved in d.c. analysis 


During use the algorithm considers all combinations of 
the nodes and the PMT; calculated is used to select the 
best set of access points. Figure 8 shows the percentage 
of unique fault signatures achieved and the effect of 
choosing the best or worst node combinations. Graphs 
of Fig. 8(a) and (b) follow similar shapes. In both cases, 
the percentage of unique fault signatures levelled off at 5 
nodes for the best combination found during the search 
algorithm. Nodes (2, 3, 4, 6 and 7) for the pre-amplifier 
and nodes (2, 3, 4, 5 and 6) for the control amplifier from 
Fig. 2, are chosen as the best combinations. A reduction 
of 4 nodes is achieved for the pre-amplifier while the 
number of nodes required for fault diagnosis of the 
control amplifier is reduced by one only from the original 
6 nodes. One of the advantages of the approach for node 
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selection adopted here is that nodes can be chosen for a 
specified percentage of unique fault signatures. For 
example, if slightly in excess of 80% of unique fault 
signatures is considered satisfactory for the ATPG of the 
preamplifier, 4 nodes are sufficient. 


15 Frequency Selection 

The application of the a.c. approach on the UUT 
involves test frequency optimization and node selection. 
The a.c. model of the circuit is shown in Fig. 4 with the 
nodes numbered sequentially. It is essential to consider 
the pre-amplifier and control amplifier together. 

The nominal frequency response of the pre-amplifier 
and control amplifier is reasonably flat over the specified 
frequency range of 20 Hz to 25 kHz. The large amount of 
computer time needed to calculate the response at many 
frequencies restricts the number which can be selected 
even before the optimization. Fight frequencies are 
chosen to demonstrate the optimization algorithm and 
are represented by the vector 


f = (20, 100, 300 Hz, 1,3,5, 10,25 kHz]" (14) 


16 An Example of Setting Up the Recognition 
Matrix 

In order to illustrate the formation of the recognition 
matrix and subsequent fault diagnosis, components R10, 
R13, R23, R25, R30, R31, R32, C15, C16, C17, C18, C19 
and C20 are varied by +50% while the other 
components are held at their nominal values. The 
deviation of the then faulty system response from the 
nominal response is evaluated at all nodes for the above 
eight frequencies. Only the amplitude ratio information 
has been used here, although phase information can also 
be added to enhance diagnosability. 

The deviation vectors evaluated are grouped to form 
the recognition matrix. Similar matrices are formed for 
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Total no. of signatures 108 
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(see fig.&) 


all nodes. The optimization algorithm shown in Fig. 7 is 
applied to each matrix and the final optimum set 


Sone [300 Hz,.1000 Hz, I0-kHz]' 


selected according to equation (13). 


(15) 


17 Node Minimization—A.C. Approach 

The minimization of nodes is carried out in the same way 
as the d.c. approach outlined in Section 8. Figure 9 
shows the percentage of unique fault signatures for the 
best and worst combinations of nodes. Nodes 12, 26, 32 
and 35 with three frequencies of equation (15) are 
selected for the application of fault diagnosis giving a 
total confidence level of 92%. 

It can be seen from Fig. 9(c) corresponding to the best 
combination of access points and test frequencies, that 
an increase of the number of access points from 4 to 8 did 
not give any noticeable increase in the percentage of 
unique fault signatures, thus obeying a typical ‘Pareto’ 
law, i.e. a large increase in resources does not result in a 
corresponding increase in effectiveness. Also, the 
asymptotic level in the percentage of unique fault 
signatures may well be achieved only for the best 
combination of access points and test frequencies. 


18 Fault Diagnosis Results 
To show that ‘the selected 4 node, 3 frequency test is 
satisfactory, 7 components are considered faulty in turn 
and a pattern recognition approach is used to achieve 
diagnosis. We proceed as follows; if 

y = Diy- -Yml (16) 


is a measurement vector due to an unknown fault and 


dik 
i=1,2,...M measurements (17) 
k — 1,2,...L faults 
100 a$ 3| S| s| 


e wn o 
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32 
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Percentage of unique fault signatures 
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(b) Test frequencies 300 Hz and | kHz 


No. of nodes used 
(c) Test frequencies 300 Hz, 1 kHz 
and 10 kHz 


Fig. 9. Percentage of unique fault signatures obtained for best and worst combinations of access points 
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is a matrix formed by grouping measurement vectors of 
known faults, then ‘distances’ between the unknown 
fault and the known faults are formulated as 


M 4 
4-|Y ba? |, k= a E (18) 
i=t 


Using the nearest neighbour rule the most likely fault is 
that with minimum d,. 

The verification of the procedure is carried out by 
imposing typical faults on the actual circuit of the pre- 
and control amplifier. In the d.c. approach faults were 
imposed on R17, R11, R20 and R18 for the pre-amplifier 
and R37, R36 and R38 for the control amplifier. The 
coded fault dictionary at the selected nodes is a priori 
formed and nearest neighbour distances calculated for 
the coded measurement vector due to the imposed fault. 
All faults imposed on the pre-amplifier and R37 and R38 
on the control amplifier were correctly classified based 
on the minimum distance. However, the 'distance' 
calculated for the measurement vector due to the 
imposed fault on R36 was found to be the same for 
representative vectors in the dictionary for R36 and R37. 

A similar procedure was employed for the verification 
of the a.c. approach and the results are tabulated in 
Table 2. Fault dictionaries for the four selected nodes 
and at each of the three chosen frequencies are formed 
and coded. If d, , d,, and d,, are the nearest neighbour 
distances for faults at frequencies f, , f; and f, respectively 
then 

d, = (dp,  d,, t d.) for k = 1, 2,. . . L fault cases 
(19) 
The most likely fault is the one for which d, is a 


minimum. 


19 Effect of Measurement Noise on Diagnosability 
In order to study the effect of measurement noise on 
diagnosability, random noise of varying levels is imposed 
on node voltage measurements and this fault isolation 


Table 2 


Some fault diagnosis results using a.c. test measurements 
on actual hardware 


Faults Diagnosis 
Imposed Obtained 
R30 correct 
C18 within groupt 
{C18; R31; R32} 
C19 correct 
R25 correct 
C20 correct 
R10 within groupt 
{R13; R10} 
R13 within groupt 
{R13; R10} 
R32 correct 


f Equal minimum distance to each fault case 
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technique applied. Random noise with nominal node 
voltage as mean, with standard deviation (ø )equal to 1%, 
2% and 3% of nominal is imposed. One hundred 
simulations were carried out for each fault imposed at 
each level of noise and the diagnosability (i.e. percentage 
of correct diagnosis) calculated. For example, a fault in 
R30 is simulated and the voltages at the selected nodes 
are calculated. Measurement noise is imposed on these 
voltages by sampling measurement noise of standard 
deviation o from a normal population. Fault diagnosis is 
now carried out. This is repeated 100 times (note that the 
imposed fault on R30 is the same) and diagnosability 
calculated as a percentage, which is (correct 
diagnosis/no. of simulations) x 100, as listed in Table 3. 


Table 3 
Effect of measurement noise on diagnosability via 
simulation 
Faults Standard Deviation 
Imposed of Measurement Noise a 
1% 2% 3% 
R30 80% 735 58% 
R15 100% 98% 97% 
R13 100% 85% 67% 
C20 100% 88% 75% 
R10 100% 83% 78% 
R31 97% 89% 70% 
C18 76% 69% 49% 
C19 75%, 66% 46% 
R32 97% 85% 87% 
C17 98% 78% 69% 
C16 100% 69% 53% 
R23 100% 80% 79% 


20 Ambiguity Groups 

In any large analogue circuit there can be many 
components, the effect of whose variations on 
performance is indistinguishable from each other. Prior 
to setting up a fault diagnosis procedure, it is useful to 
identify and group such components to form ambiguity 
groups which could then be used to enhance 
diagnosability. Complex grouping criteria have been 
discussed in the literature.?? Application of these 
techniques to the present UUT is under further 
investigation. 


21 Automation of On-line Computer 
Implementation of Fault Diagnosis Scheme 

The fault isolation procedure for the pre-amplifier and 
control amplifier has been implemented at an automatic 
test station using a minicomputer. This necessitated the 
design of special interface units.?? A block diagram of the 
whole procedure is shown in Fig.10. The computer 
program for access point selection, test frequency 
selection, and fault isolation procedure were written in 
FORTRAN. Additional to the main program, a sub- 
routine was written in MACRO?’ to read the measured 
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Fig. 10. Minicomputer-based automatic fault diagnosis facility 


analogue and digital data from the laboratory peripheral 
system to the main program. Automatic fault location is 
implemented in two phases, namely d.c. and a.c. 
approaches. If the first phase of the testing does not 
identify the fault, then the second phase of the test is 
implemented. 

The interface is a fundamental requirement of a fully 
automatic fault location scheme. Node voltages of the 
circuit are measured in sequence using software 
controlled digital outputs of the computer addressing an 
analogue multiplexer. In addition the gain setting of the 
output voltage from the interface is monitored 
automatically. 


22 Conclusions 

A fault diagnosis scheme for large-scale analogue circuits 
has been developed, based on a limited number of node 
voltage measurements to isolate drift faults to 
component level. Procedures for the selection of the best 
set of test points for d.c. and a.c. approaches have been 
developed and an optimization algorithm for the 
selection of test frequencies has been successfully applied. 
The techniques developed were applied to a hardware 
amplifier and control amplifier circuit and the whole 
procedure was automated using a minicomputer. 
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SUMMARY 


Problems inherent in conventional p-n junction solar cells 
are discussed. The principle of operation of min m.i.s. solar 
cells is then described. Among other advantages it is 
shown that the min m.i.s. solar cells are capable of giving 
a much higher Vos compared to the p-n junction solar 
cells. Problems associated with the min m.i.s. solar cells 
are also discussed. These include the incompatibility of 
light transmission and conductivity properties of the metal 
contact in the case of the contact covering the whole of 
the top surface, the contact linewidth/space ratio in the 
case of grating cells, resistivity dependence of Voc and 
degradation of cell performance in time. 
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1 Introduction 

Various methods exist for the conversion of solar energy 
to electricity, but the one that holds the best promise is 
the photovoltaic cell, usually known as the solar cell. The 
solar cell is a semiconductor junction device which 
converts sunlight directly to electricity. The 
semiconductor mostly used is single-crystal silicon, and 
the type of solar cell that has so far received the most 
attention is the p-n junction. Other semiconductors 
currently studied for solar cell fabrication include GaAs 
and CdS. Work is also in progress on polycrystalline and 
amorphous solar cells. Recently, very high efficiency 
(18:395) has been obtained for the grating type minority 
carrier  metal-insulator-semiconductor (min m.i.s.) 
solar cell. This paper discusses the principle of operation 
and properties of min m.i.s. solar cells as contenders to 
the conventional p-n-junction solar cells. 


2 p-n-Junction Solar Cell 
Figure 1(a) shows a p-n-junction solar cell. A metal grid 
on the top surface helps lower the series resistance and 
forms one of the two output contacts; the entire 
metallized bottom surface forms the other contact. This 
cell measures 2cm x2cm and has an open-circuit 
voltage of 0-56 V and a short-circuit current of about 
120 mA under AMIf conditions of illumination. Figure 
1(b) shows a cross-sectional diagram of the cell. The top 
n-diffused layer is very thin, typically 0:3 pum and 
phosphorus-doped, so that the junction is formed very 
close to the surface. Since the generation rate of electron- 
hole pairs falls off exponentially as a function of distance 
from the surface due to absorption and the effective 
carrier collection is about a diffusion length on either 
side of the junction depletion layer, a junction close to 
the surface will ensure maximum efficiency of operation. 
However, making the junction shallow increases the 
series resistance of the cell since the current generated 
has to flow to the grid contact through the n-type 
diffused layer. In order to reduce this series resistance an 
alternative to having a relatively deep junction is to 
increase the doping of the top layer. This, unfortunately, 
increases the recombination rate and thus lowers the 
carrier lifetime in this very important top region! and 
creates what is normally known as a ‘dead layer’. As this 
is a region of maximum carrier generation, the existence 
of the dead layer greatly lowers the carrier collection 
efficiency and reduces the short-circuit current /,,. The 
large recombination rate in the top layer also lowers the 
open-circuit voltage, V, through the increase in the 
saturation current, lọ. Another factor that causes an 
increase in saturation current is the narrowing of the 
semiconductor bandgap due to excessive doping. 
Another problem associated with p-n junction solar 
cells is the preferential diffusion of dopants through grain 


+ AMI: Sea level at noon on the equator. In practice the conditions are 
established by reference to a secondary standard. 
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(a) A p-n-junction solar cell 
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(b) Cross-sectional diagram of the solar cell 
Fig. 1. 


boundaries when polycrystalline materials are used. This 
preferential diffusion has the tendency of short-circuiting 
the junction formed and results in a low Ve- 

It is thus clear that there are problems inherent in the 
p-n-junction structure that are not easily resolved. The 
m.i.s. solar cell, to be discussed in the following Section, 
can in principle overcome some of the problems 
mentioned above and provide additional advantages 
such as simplicity of fabrication. 


3 Min M.I.S. Solar Cell 

In the m.is. solar cell I is a thin insulator of 
SiO, ( « 20 À). The development of the m.i.s. solar cell 
originated from the Schottky-barrier solar cell.” 

A Schottky barrier is a metal-semiconductor contact 
with rectifying properties. Because the high-temperature 
diffusion step is not required, the process is particularly 
suited to polycrystalline solar cell fabrication which has 
proved to be difficult in conventional p-n junctions due 
to junction shorting by preferential diffusion of 
impurities through grain boundaries. Figure 2 shows the 
I-V characteristics of a Schottky barrier in the dark and 
under illumination. Ideally the illuminated 
characteristics are those in the dark displaced 
downwards. In the same diagram the characteristics of a 
p-n junction in the dark as well as under illumination are 
also shown. In comparing the two sets of characteristics 
it may be seen that the Schottky barrier gives a lower I 
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Fig. 2. IV characteristics of a Schottky barrier and a p-n-junction 
solar cell in the dark and under illumination. 


and a considerably smaller V,.. The lower I. is due 
mainly to the absorption of light by the semitransparent 
metal layer which forms the Schottky barrier (Fig. 3), 
and to surface states at the metal-semiconductor 
interface. The low V,, is attributed to the basic operation 
principle of the Schottky barrier as a majority-carrier 
device. 


a Metal contact (Semi- 
transparent ) 


— —Semiconductor 


7 Metal contact (ohmic) 


Fig. 3. A Schottky-barrier solar cell. 


The incorporation of a thin insulating layer between 
the metal and semiconductor seems to have multiple 
effects towards the improvement of V," First, as shown 
in Fig. 4 the insulating layer introduces an extra voltage 
drop A@ and therefore reduces lọ. The open-circuit 
voltage is given by 


(1) 


Thus, reducing I, increases V. 

Secondly, the existence of the insulating layer 
markedly reduces the forward-bias dark current of the 
cell. This is illustrated in Fig. 5. The forward-bias current 
is mainly a majority current whose transport through the 
junction is strongly affected by the transport properties 
of the insulating layer. There is a small minority current 
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Fig. 4. Band diagram of an m.i.s. cell. 


component, however, that is controlled mainly by 
diffusion. Until the oxide is over a certain thickness, 
about 20À, the diffusion mechanism remains the 
controlling factor for the transport of minority carriers ; 
above this critical thickness tunnelling becomes an 
important controlling factor. Hence the introduction of 
the insulator can substantially reduce the majority- 
carrier component of the current and if the thickness of 
the insulator is less than, say, 20A the minority-carrier 
current is relatively unaffected. 


OÁSÀ 15À — ?22À 


Increasing 
insulator thickness 


v 


Forward - biased 


Fig. 5. Dark forward-biased IV characteristics of an m.i.s. cell as a 
function of insulator thickness. 


Thirdly, it is well known that a Schottky barrier often 
measures the same barrier height irrespective of the 
metal work function. This is because the Fermi level is 
pinned (stabilized) to a fixed level at the semiconductor 
surface by the high density of interface states between the 
metal and semiconductor.? It is also well known that an 
oxide layer provides anchorage for some of the dangling 
bonds and thus effectively frees the Fermi level which can 
then take up a position as determined by the relative 
magnitudes of the metal and semiconductor work 
functions. 

Finally, the Fermi level having been freed from the 
pinning effect of the interface states, a metal with an 
appropriate work function may be chosen so as to create 
an inversion at the semiconductor surface. This further 
reduces the majority-carrier component of the current 
and converts a principally majority-current Schottky 
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7— Inversion charge 
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Fig. 6. An m.i.s. solar cell. 
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device into a minority-current (min) m.i.s. device. The 
junction in this device is an induced one, as shown in 
Fig. 6. Since it is now a minority-current device it 
behaves in almost every respect as a conventional p-n 
junction, and as such it is expected to give the same 
efficiency of energy conversion. 

Figure 4 shows the band diagram of such a junction 
for a p-type semiconductor. A little consideration will 
show that for maximum inversion of the semiconductor 
surface a low work-function metal must be used with p- 
type semiconductors and a high work-function metal 
with n-type semiconductors. Candidates for metals in p- 
Si m.i.s. are Mg, Al and Ti, and for metals in n-Si m.i.s., 
Au, Ni and Pd. 


3.1 The Metal Contact 

In principle an m.i.s. solar cell should have the contact 
metal covering the entire surface of the cell. To let 
sufficient light through to the semiconductor the metal 
must be made very thin (< 100 A). However, such a thin 
metal film will have a high sheet resistance, and to make 
a metal film with good conductivity and reasonable 
transmission properties would require very critical 
control of its thickness and purity of the metal. At a 
thickness of less than 100:A the reproducibility can be a 
problem. Furthermore, unless well passivated, it does 
not take long for the metal to become oxidized right 
through the film. 


3.2 Transparent-oxide Heterojunction 

An alternative to the use of a metal top contact is a 
transparent oxide such as SnO,.° The electron affinity of 
SnO, is about 4-85 eV and should be sufficient to invert 
the surface of n-type silicon as shown in Fig.7. 


| | — — — —Vacuum 


Level 


Fig. 7. Band diagram of an SnO,-Si heterojunction. 


Experience has shown, however, that SnO, alone does 
not give sufficiently low sheet resistance and to improve 
the situation the film is often doped with indium. 
Unfortunately, like metal films, conductivity and 
transmission are incompatible properties, ie. a high 
conductivity invariably means low transmission and 
vice-versa, and in practice a compromise has to be 
accepted. The SnO, film, with a refractive index of 1-9- 
2:2, is itself an ideal antireflection layer between air and 
silicon and thus represents the saving of one step in the 
fabrication process. 
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3.3 Grating M.I.S. Cell 

An alternative to having the metal over the entire surface 
is the grating structure shown in Fig. 8. Underneath the 
metal grating inversion results as a consequence of the 
metal work function. The metal grating has a thickness 
of 1-4 um and thus will absorb all the light falling upon 
it. Its function here, therefore, only serves to create an 
induced p-n junction for carrier collection. The active 
part of the cell where carriers are actually generated is 
the area between the grating lines. Means must be 


— —-Mq grating 


inversion 
layer 


Substrate 


Fig. 8. An m.i.s. grating structure. 


provided, therefore, for the semiconductor in this area to 
be inverted. A convenient way to do this is to select an 
insulator to cover the semiconductor surface; the 
insulator chosen must have the right refractive index so 
as to serve as an antireflection layer and contain the right 
kind of charge so as to invert the semiconductor surface 
underneath.!9-! ! SiO seems to serve the purpose well. It 
can be evaporated in vacuum and if the rate of 
evaporation is low the resulting SiO layer is found to 
contain sufficient positive charge which is required for 
the semiconductor surface inversion. SiO is by no means 
the only insulator suitable for this application, but 
because of its refractive index of 1-9 which is ideal as an 
antireflection layer of the 4/4 type, and the fact that it can 
easily be deposited by evaporation it becomes the 
immediate candidate for such a purpose. Other 
insulators, for example, SiO;, may serve just as well. 

It may be mentioned, in passing, that a grating 
metallization on top of the transparent oxide SnO, will 
also help to reduce the series resistance of the 
heterojunction cell. 


3.4 Grating Size 

Assuming that the surface is inverted in the area 
between grating lines, the field that exists in the space- 
charge region will direct electrons generated to flow 
perpendicularly to the inversion layer. Once they are in 
the inversion layer they flow laterally to the collecting 
contact lines of the metal grating as indicated in Fig. 9. 
Power loss occurs mainly as a result of resistive effects in 
the inversion layer. Recombination with holes and 
trapping by surface states can also occur in the electron 
journey to the grating lines particularly if inversion of the 
surface is insufficient. Note that the lateral flow of 
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Fig. 9. Electron flow paths to contact lines. 


electrons increases linearly from the centre towards the 
grating lines. 

Most published results on m.i.s. solar cells give active- 
area efficiency. This is the efficiency of the solar cell 
which would be obtained in the absence of contact lines 
which only act as current collecting electrodes. These 
lines block the transmission of sunlight and are a source 
of power loss. The wider these contact lines are in 
relation to the spacing between them the lower is the 
ratio of generated power to available solar power. On 
the other hand increasing the spacing increases the 
resistive loss. A simple calculation shows that for a 
minimum power loss the spacing should be!? 


m (my 
Psd mp 

where W is the width of the contact lines, p, is the sheet 
resistivity of the inversion layer in O/square, and V4, and 
Jmp are the generated voltage and current density under 
maximum power conditions. If W= 10pm, 
p, = 5x 10* Q/square, Vap = 450 mV and 
Jap = 30 mA/cm?, 5-121 um, which is about that 
presently used for m.i.s. grating cells (order of 100 uim). 


Q) 


3.5 Antireflection Coating 

Most semiconductors have a refractive index greater 
than 3, and therefore a high reflectance. An antireflection 
layer is usually formed by depositing over the 
semiconductor surface a coating of material which is 
transparent to sunlight in the spectrum of interest and 
which has a refractive index n, such that ng « n, « n. It 
may be shown that the reflectance is a minimum if 


(3) 
Equation (3) is the condition for 4/4 matching, and the 
thickness t, of the antireflection layer is given by 
nyt, = A/4 (4) 

To match silicon (n = 3-6) to air (ng = 1), n, should be 
about 1-9. Equation (4) gives a thickness of about 79 nm, 
taking 600 nm as the midpoint of the useful spectrum. 
The reflectance over the spectrum of such a system is 
shown in Fig. 10.'3 

Better results can be obtained if two antireflection 
coatings are used. The result is also shown in Fig. 10 for 
a particular combination of two antireflection layers. 


nl = non 
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Fig. 10. Reduction of reflectance by antireflection layers. 


For use in the m.i.s. grating solar cell the antireflection 
coating has the additional responsibility of inducing 
inversion on the semiconductor surface. SiO has been 
found to be suitable for p-type material, i.e. it will induce 
an n-type surface, while Al,O, is suitable for n-type 
material, i.e. it will induce a p-type surface. SiO can be 
deposited by evaporation in vacuum while Al,O, can be 
formed by anodization of aluminium deposited also by 
evaporation in vacuum. 


3.6 Textured Cells 

Another way to reduce the reflectance is by texturing the 
semiconductor surface.|*!5 This is done by a 
preferential etch which etches silicon more rapidly in one 
crystallographic direction. As a result pyramids are 
formed and light undergoes multiple reflections as shown 
in Fig. 11. Together with an antireflection coating the 
loss of light due to reflection can be kept below 3%. 
Furthermore, light enters the surface at an angle and is 
therefore absorbed close to the surface. This should 
improve the carrier collection efficiency, particularly for 
light of longer wavelengths. This is an advantage 
especially for the m.i.s. solar cell since the *dead layer' as 
found in the conventional p-n junction cell has been 
eliminated. 


Incidence 


Fig. 11. Transmission and reflection of light at a textured surface. 


The texturing, however, means an additional 
processing step which requires careful handling and thus 
increases the processing cost. Furthermore, because of 
the textured surface contact grating lines become longer 
by a factor of X 3. This increases the series resistance of 
the cell unless the grating lines are made correspondingly 
thicker or wider. 
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3.7 Effects of Substrate Resistivity on Voc 

Figure 12 shows the highest V,, so far attained.!$ The 
considerably higher V,. obtained in miis. cells is a 
consequence of the elimination of the ‘dead layer’ that 
exists in p-n junction solar cells. 


o00MIS 
aœ o pn junction 


105 21-5 -2105-020 
Resistivity (.^.- cm) 


Fig. 12. V,, as a function of substrate resistivity. 


In general, since the saturation reverse current I is 
given by 
q, — GADani 
qc eu 
LN, 


where A is the area of the cell, D, the diffusion constant 
for electrons, n; the intrinsic concentration, N, the 
substrate doping concentration and L, the electron 
diffusion length, increasing N, decreases Jo. Therefore, 
as seen from equation (1), VW. increases with N4. 
However, the lifetime of minority electrons falls rapidly if 
N, is increased beyond a certain limit, and J will 
increase. Figure 12 shows that maximum V,, occurs at a 
substrate resistivity of about 0-1 Q cm, corresponding to 
N, > 5x TOP cme, 

Another factor that limits V,. in min m.i.s. solar cells is 
the degree of inversion that may be obtained with a given 
substrate resistivity. The inversion charge Q, is given 
by?’ 


(5) 


Qn = Q.. — (4qe £o N A Dg, (6) 
where Q, is the fixed positive oxide charge responsible 
for the inversion, and $, is the Fermi potential given by 

kT Na 


Ep ee 
Ore q i ni 


(7) 


Therefore, the higher N, is, the lower is Q, and hence /,,, 
the short-circuit current. Again, as seen from 
equation (1), a decrease in L, will also result in a 
decrease in V, In addition, a high N, means a narrow 
depletion width, leading to a poor response for the long 
wavelengths, unless a back surface field (b.s.f.) is 
incorporated. 


3.8 Efficiency n 

As explained in Section 2, the main problem with the 
conventional p-n junction solar cell is the existence of the 
‘dead layer’ which is associated with the heavy doping of 
the front emitter region. The min m.i.s. solar cell, on the 
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other hand, eliminates this heavily-doped layer. Instead, 
an inversion layer exists immediately at the surface of the 
semiconductor. The induced junction is electronically 
identical to the ideal p-n junction and thus should have 
similar properties.S However, compared with the 
practical p-n junction the elimination of the ‘dead layer’ 
in the min m.i.s. solar cell greatly enhances the carrier 
collection efficiency leading to an increase in I,,. It also 
lowers the saturation current, Ig, leading to an increase 
in V4. Since the induced junction is actually at the 
surface, there is a marked improvement in the short- 
wavelength response as compared to a p-n-junction 
cell! (approximately 85% increase in response is 
obtained at 370 nm). 

Another factor that helps increase the efficiency of the 
min m.i.s. solar cell is the elimination of the diffusion 
step in the fabrication process. This saves the cell from 
being subjected to high-temperature treatment which is 
known to cause deterioration in carrier lifetime.'? 

Recent experimental results! give an active-area 
efficiency of 18:3% for the min m.i.s. solar cell without 
any back surface field. This compares very favourably 
with that of the more complicated black cell'? whose 
active-area efficiency is 18:5% or the p* —n— n* b.s.f. 
cell?° whose active-area efficiency is estimated as 18:697. 
It is worth noting the very high V4, obtainable from the 
min m.i.s. solar cell. At very high V4, (> 650 mV), surface 
states become important as they cause excessive surface 
recombination. There are thus prospects for further 
increase of efficiency in the min m.i.s. cell, for example, 
by incorporating a back surface field in the base region, 
and by proper surface treatment. These effects are 
presently being investigated. 


3.9 Degradation Problems 

Figure 13 shows the characteristics of a min m.i.s. solar 
cell with magnesium contact lines immediately after 
fabrication and a week later. The contact-line spacing of 
this cell has been deliberately made very wide (about 
2:6 mm) in order to reveal specifically some degradation 
effects. It is clear that in the period of one week severe 
degradation has taken place. 


V (VOLT) 
Fig. 13. IV characteristics of a coarse-grating cell. 


A possible mechanism leading to the degradation is 
that of the m.i.s. contact-line structure. This is due to the 
reduction of the thin oxide layer by the overlying metal 
contact. High-temperature life-testing of Al/SiO,/p-Si 
solar cells?! has shown that below 200°C this 
mechanism proceeds at a negligible rate. Mg/SiO,/p-Si 
solar cells are presently under investigation and 
preliminary results seem to suggest that the problem is 
not serious. 

Another possible mechanism for degradation is caused 
by the ageing effect of the insulator-semiconductor 
interface leading to an increase in surface-state density. 


TWO MONTHS LATER 


0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 


E - E, (eV) 
Fig. 14. N,, curves of an Al-SiO-Si structure. 


Measurements on some m.i.s. capacitors with SiO as the 
insulator have shown that in a period of two months an 
order of magnitude increase in surface-state density is 
possible (Fig. 14). However, if the semiconductor surface 
is inverted most of the surface states will be occupied 
under steady-state conditions and unable to cause 
degradation on a magnitude as shown in Fig. 13. 

The most likely mechanism leading to degradation 
originates from the weakening of the inversion layer at 
the semiconductor surface underneath the SiO in 
between contact lines. This inversion layer assists the 
collection of minority electrons in at least three ways. 
First, the inversion layer provides a path of low 
resistance for the lateral movement of electrons. 
Secondly, immediately below the inversion layer is the 
depletion region whose width is a maximum at strong 
inversion with the associated large electric field. À large 
electric field is essential for the separation of the electrons 
from the holes and for driving the former towards the 
semiconductor surface. Thirdly, electrons from the 
inversion layer help immunize the effects of surface states 
present at the insulator-semiconductor interface, which 
otherwise would capture the  optically-generated 
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Fig. 15. H.f. CV characteristics of an M-SiO-Si capacitance. 


minority electrons. If for any reason inversion disappears 
altogether, leaving behind only the depletion layer, 
electrons arriving at the semiconductor surface under the 
depletion-layer field will suffer serious resistive loss in 
their lateral movement to the contact lines due to 
recombination at the surface states as well as low 
electron concentration at the surface. The collection of 
minority electrons may then be limited to regions within 
a diffusion length on either side of the contact lines. 
Evidence of the disappearance of the inversion layer on 
either side of a contact line is clearly shown in Fig. 15 
where CV curves of an m.i.s. capacitor are shown with 
time elapsed after fabrication. The insulator of the m.i.s. 
capacitor is SiO, of about 800A in thickness. The 
frequency of measurement is 1 MHz. Low-frequency CV 
characteristics are evidence of the existence of an 
inversion layer in regions of the semiconductor 
surrounding the metal contact, the so-called a.c. 
conductance effect.?? It can be seen that in about two 
days after fabrication the CV characteristics have 
changed from the low-frequency type to the high- 
frequency type, indicating a considerable weakening of 
the inversion layer surrounding the contact. Unless this 
problem can be overcome line spacing should be made 
less than a diffusion length. In top-grade silicon a 
diffusion length in the region of 130-180 jum may be 
expected. This compares favourably with the line spacing 
of about 100 um presently used for m.i.s. solar cell 
contact grating. For wider spacings degradation may be 
expected in time as the inversion layer weakens. 

The weakening or disappearance of the inversion layer 
may be caused by the neutralization of the positive 
charge in the antireflection coating or the accumulation 
of negative charge on the top surface of this coating. 
Experiments have shown that the accumulation of 
negative charge is the main reason responsible for the 
degradation.?? The problem, however, may be overcome 
by appropriate passivation.?* One of the methods tried is 
to cover the cell with a layer of concentrated non-ionic 
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surfactant C0630.t This is a colourless and very viscous 
liquid. Experiments on C0630-passivated fine-grating 
cells showed that over a period of three months, I 
remained practically constant at its maximum value 
while for the unpassivated cell, /,, fell by at least 5%. 
C0630-passivated coarse-grating cells showed a variation 
in I, of about + 107; in a period of 8 months while in the 
unpassivated cell J,, fell to about 20% of its starting 
value. Other non-ionic surfactants may also give 
acceptable results. For example, C0990, a solid non-ionic 
surfactant at room temperature, may be applied to the 
surface of the cell by heating and spinning. However, 
because C0990 is an off-white wax it causes some loss of 
efficiency in photovoltic operation due to light 
absorption. It has also been found possible to use 
mixtures of C0630 and C0990 which retain good 
transmission properties while still remaining a solid at 
room temperature. Research is at present directed to 
finding an insulator or combination of insulators which 
would serve the combined purposes of antireflection, 
inversion charge inducement and passivation. It may be 
mentioned, in passing, that problems associated with 
this type of degradation will also occur in other solar 
cells which use antireflection coatings as charge inducing 
layers!? so that similar considerations will apply. 


4 Conclusions 

In conclusion it is useful to summarize the advantages 
and disadvantages of the min m.s. solar cells as 
compared to the conventional p-n-junction solar cells. 

l. The min m.i.s. solar cell is a relatively simple 
structure to fabricate and economic studies?? have 
shown lower cell processing costs. 

2. Elimination. of the heavily doped top layer 
removes the ‘dead-layer’ that exists in the conventional 
p-n-junction, thus increasing the carrier lifetime and 
collection efficiency in the surface region of the cell where 
the generation rate is a maximum, especially for the 
ultra-violet end of the spectrum. In addition, the 
potential barrier of the min m.i.s. cell is right at the 
semiconductor surface and this further improves the 
collection efficiency of carriers generated by light of short 
wavelengths. The combined effect is an increase in Ie 
and Vie- 

3. Elimination of the diffusion step in the fabrication 
process saves the cell from being subjected to high- 
temperature treatment which is known to cause 
deterioration in carrier lifetime. This should further help 
improve the efficiency of the min m.i.s. solar cell. 

4. Elimination of the diffusion step makes the 
min m.i.s. cell an ideal structure for polycrystalline 
materials since there is no problem of preferential 
diffusion of dopants through grain boundaries which 
tends to short circuit the junction. The use of 


+ Surfactant C0630: trade name Antarox, manufactured by GAF 
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polycrystalline materials should result in a reduction of 
cost. 

5. Min mi.s. solar cells are faced with the problem of 
finding a top contact material that is sufficiently 
transparent and at the same time has a high 
conductivity. The grating min m.i.s. solar cell, on the 
other hand, tends to be more susceptible to degradation. 
Three degradation mechanisms have been identified. 
The principal mechanism found to be responsible for the 
degradation is the accumulation of negative charge on 
the outer surface of the antireflection coating. However, 
recent research to overcome this problem seems very 
promising, and it may be concluded that when this 
problem is entirely solved, the grating min m.i.s. solar 
cell will be a real contender to the conventional p-n- 
junction solar cell for the direct conversion of solar 
energy into electricity. 
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SUMMARY 


The paper describes various detection processes for a 
modem that has been designed for operation at 9600 bit/s 
over the public switched telephone network. The modem 
is a synchronous serial system, using a 16-point q.a.m. 
signal, with near-maximum-likelihood detection at the 
receiver. The detection processes are modifications of a 
technique that has recently been developed for use with 
severely distorted digital signals, the modifications leading 
to much simpler implementation. Results of computer- 
simulation tests are presented, comparing the tolerances of 
the various detectors to additive white Gaussian noise, 
when the data are transmitted over different telephone 
circuits. The effects, on the tolerance to noise of a 
detector, of inaccuracies in the estimates of the level and 
carrier phase of the received signal are also studied. 
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1 Introduction 


A telephone circuit in the switched telephone network 
may have any one of a very wide range of attenuation- 
frequency characteristics and any one of a 
correspondingly wide range of group-delay frequency 
characteristics, giving an available (or usable) bandwidth 
that may lie anywhere in the range from about 1500 Hz 
to 3000Hz.! It is evident therefore that for the 
satisfactory transmission of data at a rate as high as 
9600 bits/second over the switched telephone network, 
the receiver must be adaptive in the sense that it takes 
full account of the distortion introduced into the received 
data signal. The conventional approach to this problem 
is to use an adaptive non-linear (decision-feedback ) 
equalizer that is adjusted to minimize the mean-square 
error in the equalized signal.” ? Such systems are known 
to operate satisfactorily over telephone circuits at rates of 
up to 4800 bits/second. At higher transmission rates, 
however, satisfactory operation is not always obtained 
over the poorer telephone circuits. 

An alternative approach that has recently been con- 
sidered is the use of the Viterbi-algorithm detector. ? 
This selects as the detected message the possible 
sequence of transmitted data-symbols (signal-element 
values) for which there is the minimum mean-square 
difference between the samples of the corresponding 
received data signal, for the given signal distortion but in 
the absence of noise, and the samples of the signal 
actually received. When the data signal is received in the 
presence of stationary additive white Gaussian noise, 
giving statistically independent Gaussian noise samples 
at the detector input, this is a process of maximum- - 
likelihood detection; when the transmitted data-symbols 
are statistically independent and equally likely to have 
any of their possible values, the detection process 
minimizes the probability of error in the detection of the 
received message.? It is assumed here that the received 
signal is sampled at the Nyquist rate, so that the 
corresponding samples contain all the information in this 
signal. 

Unfortunately, when the sampled impulse-response of 
the channel contains a large number of components 
(non-zero samples), which is the case in the particular 
application considered here, the  Viterbi-algorithm 
involves both an excessive amount of storage and an 
excessive number of operations per received data 
symbol. One approach for overcoming this difficulty is to 
use a linear feedforward transversal filter at the detector 
input to reduce the number of components in the channel 
sampled-impulse-response. The filter is adjusted 
adaptively to give a ‘desired’ sampled impulse-response 
for the channel and filter, which has a given small 
number of components and may or may not be fixed.2""! 
The disadvantage of this arrangement is that, for some 
telephone circuits, the linear filter may equalize some of 
the amplitude distortion introduced by the telephone 
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circuit and under these conditions maximum-likelihood 
detection is no longer achieved, leading to an inferior 
performance.? The linear filter should ideally perform 
the function of a ‘whitened matched-filter',$:! ^ which 
ensures that true maximum-likelihood detection is 
achieved by the Viterbi-algorithm detector. This filter 
contains a linear filter matched to the channel, followed 
by a transversal filter, the two filters being separated by a 
sampler. The sampler samples the output signal from the 
matched filter, once per data symbol, the sampling phase 
being appropriately adjusted to maximize the 
signal/noise ratio in the samples. The zeros (roots) of 
the = transform of the sampled impulse-response of the 
channel and matched filter occur in complex-conjugate 
reciprocal pairs, and it is assumed here that there are no 
zeros exactly on the unit circle in the - plane (that is, 
have an absolute value of unity). The samples are fed to a 
linear feedforward transversal filter which removes all 
zeros of the = transform of the channel and matched 
filler that lie outside the unit circle, leaving the 
remaining zeros unchanged. Thus the z transform of the 
channel and whitened matched filter has no zeros outside 
the unit circle. 

In the ideal arrangement, where the highest possible 
transmission rate is required over a baseband channel 
which is bandlimited and introduces stationary additive 
white Gaussian noise at its output, the real-valued 
baseband data signal is transmitted over the channel at 
the Nyquist rate. The receiver input filter is now a 
lowpass filter having a flat amplitude response and linear 
phase characteristic over the passband and a cut-off 
frequency (in Hz) at half the signal element rate. The 
output signal from the lowpass filter is sampled once per 
data symbol (at the Nyquist rate). The samples of the 
data signal contain all the information in the data-signal 
waveform and the corresponding noise samples are 
statistically independent Gaussian random variables 
with zero mean and fixed variance. 

Under these conditions the linear matched filter and 
sampler that form the first part of the whitened matched 
filter can be replaced by the given receiver input filter and 
sampler, followed by a linear feedforward transversal 
filter, which is such that the resultant sampled impulse- 
response (for the receiver filter, sampler and transversal 
filter) is the same as that of the original matched filter 
and sampler. Since the receiver input filter and sampler 
introduce no information loss (there being also no 
bandlimiting of the data signal by the receiver filter), they 
may now be considered as part of the channel, whose 
sampled impulse-response (for a given sampling phase) 
remains unchanged. The transversal filter therefore 
becomes the matched filter for this channel. The sampled 
impulse-response of the transversal filter is just the time 
reverse of that of the channel, and, in the more general 
case with complex-valued signals, it becomes the complex 
conjugate of the time reverse. The zeros of the z 
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transform of the matched filter are the complex 
conjugates of the reciprocals of the zeros of the - 
transform of the channel. Thus, for every zero of the 
channel z-transform, the matched filter adds another 
zero at the complex conjugate of the reciprocal value of 
z, such that the zeros of the resultant z transform occur in 
complex-conjugate reciprocal pairs. 

Following the transversal matched filter just described 
is the linear feedforward transversal filter previously 
mentioned, that removes all zeros of the = transform of 
the channel and matched filter that lie outside the unit 
circle. The two transversal filters operate with the same 
sampling interval and phase and can clearly be combined 
to give a single linear feedforward transversal filter. It 
can be seen that this filter must be such that, in the z 
transform of the sampled impulse-response of the 
channel and filter, all zeros of the z transform of the 
channel, that lie outside the unit circle in the z plane, are 
replaced by the complex conjugates of their reciprocals, 
all remaining zeros being left unchanged. Thus all zeros 
of the z transform of the channel and filter lie inside the 
unit circle. It has been shown that this filter is an all-pass 
network (with a flat amplitude response) that performs 
an orthogonal transformation on the received signal.? 
The filter removes all phase distortion introduced by the 
channel and then processes the signal, without changing 
the amplitude distortion introduced by the channel, in 
such a way as to place all zeros of the z transform of the 
channel and filter inside the unit circle. Thus, whenever 
the channel introduces any amplitude distortion, the 
channel and filter introduce the same amplitude 
distortion together with some phase distortion, the latter 
being determined entirely by the amplitude distortion. 

The linear transversal filter of the whitened matched- 
filter just described is identically the same as the linear 
feedforward transversal filter that forms the first part of a 
conventional non-linear (decision-feedback) equalizer, 
where this is fed with the given sampled signal and is 
adjusted to minimize the mean-square error (and hence 
maximize the signal/noise ratio) in its output signal, 
subject to the exact equalization of the channel.? It 
follows that the linear transversal filter of the whitened 
matched filter may be adaptively adjusted for the 
channel to give approximately the required response, at 
least at high signal/noise ratios (where near perfect 
equalization of the channel is achieved by the 
corresponding adaptive equalizer), using basically the 
same simple arrangement as that involved in a 
conventional adaptive non-linear equalizer.? The signal 
at the output of the linear filter is fed to the maximum- 
likelihood detector, and the feedback filter gives an 
estimate of the sampled impulse-response of the channel 
and linear filter, which is required by the detector. The 
details of the arrangement will be described later. In 
practice, the linear filter usually introduces some 
constant gain or attenuation, which however has the 
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same effect on the noise as on the data signal, and so does 
not change in any way the statistical relationship 
between the noise and data signal. Furthermore, when 
an allowance has been made for the change in level, the 
noise components at the filter output have the same 
statistical properties as those at its input, being 
therefore statistically independent and with a fixed 
variance. This means that a Viterbi-algorithm detector 
operating at the output of the filter has the same 
tolerance to the Gaussian noise as the appropriate 
Viterbi-algorithm detector operating at its input, 
although generally involving much less complex 
equipment. Unfortunately, since the linear filter does not 
now remove any of the amplitude distortion introduced 
by the channel, it does not always, or even in general, 
achieve a sufficient reduction in the number of 
components of the resultant sampled impulse-response, 
to permit the practical use of a conventional Viterbi- 
algorithm detector at its output. 

It follows that the Viterbi-algorithm detector must be 
replaced here by a near-maximum-likelihood detector 
that uses very much less storage and requires far fewer 
operations per received data-symbol, for the given 
number of components in the sampled impulse-response 
of the channel. Various systems of this type have been 
developed and tested by computer simulation, with very 
promising results.?: 12-2 

An alternative approach that has recently been studied 
is to discard the linear filter altogether, so that the near- 
maximum-likelihood detector operates directly on the 
original sampled signal, at the output of the sampler.!4 
An estimator is, of course, required here to determine the 
sampled impulse-response of the channel, but this can be 
implemented relatively simply, requiring a transversal 
filter with only as many taps as there are components in 
the sampled impulse-response of the channel.!>-!® 
However, computer-simulation tests (whose results have 
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not been published) have shown that in the presence of 
very severe phase distortion, such as can occur over the 
poorer switched telephone circuits, these systems require 
a considerable amount of storage and a large number of 
operations per received data-symbol, in order to achieve 
near-maximum-likelihood detection. The most cost- 
effective arrangement for the given application is 
therefore to use an adaptive linear filter of the type 
previously described, connected between the sampler 
and detector and acting as a whitened matched filter. 
The adaptive filter has the important property that it (at 
least approximately) removes all the phase distortion 
introduced by the channel and then adjusts the sampled 
impulse-response of the channel and filter into a form 
that is ideally suited to a near-maximum-likelihood 
detector of the type considered here. '? 

The paper describes the application of the methods 
involved in some recently developed detection processes 
to the design of a q.a.m. (quadrature amplitude 
modulated) system operating at 9600 bits/second with a 
16-point q.a.m. signal. Details are given here of some 
novel detection processes that avoid the high complexity 
involved with a 16-point signal in the previously 
described processes. Results of computer-simulation 
tests are then presented, comparing the performances of 
the various systems, when operating over different 
telephone circuits, and showing also the effects of 
inaccuracies in the estimates of the level and carrier 
phase of the received signal. 


2 Model of Data Transmission System 
The model of the data-transmission system is shown in 
Fig. 1. This is a synchronous serial l6-point q.a.m. 
system. The information to be transmitted is carried by 
the data-symbols (s;!, where 

s= tatjp (1) 


andj = y —1,a = 1 or 3 and fl. 5 1 or 3. The (sj) are 


Receiver 
filter 


| Fig. 1. Data-transmission system. 
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statistically independent and equally likely to have any 
of their 16 possible values. It is assumed that s; = 0 for 
i € 0, so that sjó(r — iT) is the ith signal element at the 
input to the transmitter filter. The transmission path is a 
linear baseband channel that includes a telephone circuit 
together with a linear q.a.m. modulator at the 
transmitter and a linear q.a.m. demodulator at the 
receiver. The basic structure of the modulator and 
demodulator has been described elsewhere.'’ The signals 
transmitted over the ‘in-phase’ channel of the q.a.m. 
system are represented by real-valued quantities, and the 
signals over the 'quadrature' channel by imaginary- 
valued quantities, to give a resultant complex-valued 
baseband signal at both the input and output of the 
transmission path in Fig.l. The carriers of the 
transmitted signals in the in-phase and quadrature 


channels are taken to be V/A cos 2nf.t and 


—,/2sin 2nf.t, respectively, where f, = 1800.'" The 
transmitter filter, transmission path and receiver filter 
together form a linear baseband channel whose impulse 
response is the complex-valued function q(t), with an 
effective duration of less than (g + 1)T seconds, where g is 
the appropriate integer. It is assumed for the purpose of 
this study that q(t) is time invariant over any one 
transmission. The various types of additive and 
multiplicative noise normally introduced by telephone 
circuits are neglected here, and it is assumed that the 
only noise is stationary white Gaussian noise with 
complex values, zero mean and a flat (frequency 
independent) power spectral density, which is added to 
the data signal at the output of the transmission path, to 
give the complex-valued Gaussian noise waveform c(t) at 
the output of the receiver filter. Although telephone 
circuits do not normally introduce significant levels of 
Gaussian noise, the relative tolerance of different data- 
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transmission systems to white Gaussian noise is a good 
measure of their relative overall tolerance to the additive 
noise actually experienced over telephone circuits.’ 

The waveform at the output of the receiver filter is the 
complex-valued signal 


l 
pt) 2 Y, sq(t—iT)+ v(t) (2) 


i=l 

where l is some large positive integer. The waveform p(t) 
is sampled once per data symbol at the time instants {iT} 
to give the received samples {p;}, which are fed to the 
adaptive linear filter (Fig. 1). This filter is taken to have 
been adjusted exactly as previously described so that it 
acts as a whitened matched filter, the only difference 
being that the filter now operates with complex-valued 
samples in place of the real-valued samples previously 
assumed. Thus the sampled impulse-response of the 
linear baseband channel and adaptive linear filter 
(Fig. 1) is given by the (g+1)-component vector 


Y= Yo Yi- -Yy (3) 


which has complex-valued components and a z 
transform with all its zeros inside or on the unit circle in 
the z plane. Furthermore, when an allowance has been 
made for the change in level, the noise components {w;} 
at the output of the filter have the same statistical 
properties as the noise components {v;} at its input, 
where v; = v(iT). The delay in transmission over the 
baseband channel and adaptive filter, other than that 
involved in the time-dispersion of the received signal, is 
for convenience neglected here, so that yọ # 0 and y; = 0 
for i < 0 and i > g. Thus the sample value at the filter 
output, at time t = iT, is the complex-valued quantity 
g 


r= p Si-nYn T Wi (4) 


h=0 


2: 
[aca dh 
CG) 


Estimator for response of 
channel and linear filter 


Fig. 2. Adaptive adjustment of linear filter and estimator. 
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where the real and imaginary parts of the [wj] are 
statistically independent Gaussian random variables 
with zero mean and variance o?. 

The adjustment of the adaptive linear filter and the 
estimation of Yis achieved by the arrangement shown in 
Fig. 2. It can be seen that this is a simple development of 
a conventional adaptive non-linear (decision-feedback ) 
equalizer for the case where there is a delay in detection 
(of two sampling intervals in this particular example). 
Clearly, when there is no delay in detection, the 
arrangement degenerates into a conventional equalizer. 
The main tap of the -tap adaptive linear feedforward 
transversal filter, is close to the last tap, for all channels. 
A square marked T is a store that holds the 
corresponding sample value, the stored values being 
shifted one place in the direction shown during each 
sampling interval. A square marked X is an accumulator 
that sums the input samples, and * indicates that the 
corresponding sample value is replaced by its complex 
conjugate. a and b are small positive constants, and @ 
radians is the carrier phase correction necessitated by a 
phase error in the coherent demodulator. œ is 
determined by the carrier phase control circuit, which is 
not shown in Figs. ! and 2. The adaptive adjustment of 
$, together with that of the transversal-filter tap-gains in 
Fig. 2, is arranged to be such that yg = l, and the 
method of adjustment of each transversal filter is the 
same as that in a conventional equalizer.? Further details 
of this and of the determination of $ are however beyond 
the scope of the paper. 

Except where otherwise stated it will be assumed that 
the adaptive linear filter is correctly adjusted and that the 
estimate of — Y, | which is the sequence 
Yo. Vis --- y, Stored in the channel estimator, is also 
correct. The estimate of Y, together with a prior 
knowledge of the possible values of s;, are used by the 
near-maximum-likelihood detector and intersymbol- 
interference canceller. 

The detector (Fig. 1) operates on its input samples (rj) 
to give the finally detected — data-symbols 
Isi). s; being determined after the receipt of r/,,,, where 
n « g, so that there is a delay in detection of n sampling 
intervals. The — intersymbol-interference canceller 
removes from the samples {r;} estimates (detected values) 
of all components involving data symbols [sj] whose 
final detected values {s;} have already been determined. 
Thus the intersymbol interference canceller operates on 
r; to give 


g 
i 
Fi =r 


i y SI YA (5) 


h-ntl 


When si-n = s; 4. for all {h}, as will be assumed for the 
present, 


r= Y seb; (6) 
cO 
so that the sampled impulse-response of the channel and 
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adaptive filter has been reduced from g--1 to n+] 
components, thus greatly simplifying the detection 
process when n « g. 

Let S,, Ri, and W, be the k-component row vectors 
whose ith components are s;, r; and wj, respectively, for 
i=1,2,...,k. Also let X,, Z, and U, be the k- 
component row-vectors whose ith components are x, z; 
and u;, respectively, for i= 1,2, .. ., k, where x; has one 
of the 16 possible values of s;, 


zi = Xi khYh (7) 
h=0 


and u; is the possible value of w; satisfying 
rj — zu, (8) 


In the k-dimensional complex vector space containing 
the vectors R}, Z, and U,, the square of the ‘unitary’ 
distance between the vectors Rj, and Z, is 


IU, = tu? +u? + LL. + dul? (9) 


where |u;| is the absolute value (modulus) of u;. 

When all real and imaginary components of the {w;} 
are statistically independent and with a fixed variance, 
the maximum-likelihood vector X, is its possible value 
such that |U,|? is minimized. Under the assumed 
conditions this X, is the possible value of S, most likely 
to be correct. 


3 System A 

The detection process used here is a modification of an 
arrangement previously!? referred to as System 2 and it 
operates as follows. Just prior to the receipt of the sample 
r, at the detector input, the detector holds in store m n- 
component vectors {Q,_,}, where m is a multiple of 4, 


Oii x e m1 Xie (10) 


and x; is as previously defined. Each vector Q,_, is 
associated with the corresponding cost |U, _ ,|? (eqn. (9)), 
in the evaluation of which it is assumed that s; = s; for 
every i. This implies that the intersymbol interference 
canceller operates in the ideal manner and means, of 
course, that there is an inaccuracy in |U,_,|? whenever 
one or more of the {s;! are incorrect. 

On the receipt of rą, each of the stored vectors (Q, ,] 
is expanded into four (n 4-1) — component vectors {P,}, 
where 


(11) 


The first n components of P, are as in the original vector 
Q,-, and the last component x, has the four different 
values +2+ 2j. The cost associated with each vector P, 
is evaluated as 


P= Xp hae >? ME 


cy = |U, 4 Jr. — 25 X, nal? (12) 


h=0 
For each of the four possible values of x,, the detector 
then selects the 1m vectors {P,} having the smallest costs 
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(c.l, to give a total of m selected vectors together with 
their associated costs. Each selected vector P, is next 
expanded into four vectors {P,} where the first n 
components are again as in the original vector Q, .. ,, and 
to the given value of the last component x, (now 
«2-4 2j) are added the four different values + 1+). The 
cost associated with each expanded vector is evaluated as 
n 
IU,? = [U, ^ Ir — > Xy - aol? (13) 
b=0 
The detected data-symbol s;,_, is then taken as the value 
of x,_, in the vector P, with the smallest cost, and the 
first symbol x,_, in each vector P, is discarded to give 
the corresponding vector Q,, which is, of course, 
associated with the same cost |U,|?. Finally, m vectors 
{Q,} are selected from the 4m vectors (Q,j. as follows. 
When m = 32, the detector selects, for each of the 16 
possible values of x,, the two vectors {Q,} that have the 
smallest costs and differ in the value of x,_,. When 
m = 16, the detector selects, for each of the 16 possible 
values of x,, the vector Q, having the smallest cost. The 
selected vectors are now stored together with their 
associated costs. These arrangements can be 
implemented simply and ensure that all selected vectors 
are different, thus preventing any ‘merging’ of the stored 
vectors.'? 
The technique just described is an arrangement of 
‘double expansion’ in which it is assumed that 


Xk = Xak t Xk (14) 


where 

(15) 
and 

(16) 


so that x, is treated as the sum of two separate 4-level 
data-symbols x, , and x, ,. In the expansion of the m 
stored vectors {Q,_;}, x, , is set to zero so that x, is 
treated as though it were x, ,. In the expansion of the m 
selected vectors {P,}, x, is taken as Xa, x+ x, ,, the value 
of x, , for each vector being, of course, that determined 
in the first process of expansion and selection. 

The arrangement relies on the fact that if some 
complex number q, is at a smaller unitary distance from 
a given one of the 4 possible values of x, , than from the 
3 remaining possible values of x, ,, then it is also at a 
smaller unitary distance from x, , + Xp, x for the given 
x, , and for one of the 4 possible values of x, ,, than it is 
from the remaining 15 possible values of x, , + x, ,. This 
implies that the possible value of x, closest to y, may be 
determined in two successive operations: first the 
selection of the possible value of x, , closest to q,, and 
then, for the selected value of x, ,, the determination of 
the possible value of x, , such that x, , + xy, , is closest to 
gi 
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4 System B 

This is a development of an arrangement previously!? 
referred to as System 1 and it operates as follows. Just 
prior to the receipt of the sample r} at the detector input, 
the detector holds in store m n-component vectors 
(Q,.,! together with the associated costs (|U,. |^] 
(eqns. (10) and (13)). On the receipt of r, the detector 
expands each vector Q,. , into the corresponding 16 
vectors {P,} (eqn. (11)) having the 16 possible values of 
x,(given by ta+jf, fora = 1 or3and f = 1 or3) and it 
evaluates the cost |U,|? for each of the 16 vectors {P,;. It 
then selects the vector P, with the smallest cost and takes 
the detected data-symbol s;,_, to have the value of x, _, 
in the selected vector P,. All vectors {P,} for which 
Xk-n É 5,5, are now discarded, and the first component 
of each of the remaining vectors {P,} is omitted to give 
the corresponding n-component vectors {Q,}. The one of 
these vectors derived from the vector P, with the smallest 
cost |U,|? is the first selected vector Q,. The detector then 
selects from the remaining vectors {Q,} the m— 1 vectors 
associated with the smallest costs, to give a total of m 
vectors {Q,} and their associated costs, which are stored. 
The discarding of the given vectors {P,} prevents the 
merging of the stored vectors, since it ensures that, if 
these are all different at the start of transmission, no two 
or more of them can subsequently become the same. 


5 System C 

The technique of double expansion, applied in System A, 
is here applied to System B, to give an arrangement 
involving less storage and a smaller number of 
operations in a detection process than does System B. 
Following the expansion of the m stored vectors {Q,_,} 
into 4m vectors {P,!, for which x, = +2+42j, and the 
evaluation of the associated costs {c,}, as in System A, 
the detector selects the m vectors {P,; with the smallest 
costs, regardless of the values of any of their {x,. 
Following the expansion of the latter vectors into 4m 
vectors {P,' and the evaluation of their costs {|U,|7}, 
again as in System A, the detector selects the vector P, 
with the smallest cost and takes the detected data- 
symbol s,_, to have the value of x,., in the selected 
vector P,. All vectors {P,} for which x,_, Æ Sk-n are 
now discarded, and the detection process proceeds as in 
System B. 


6 System D 

This is an alternative approach to that in System C 
towards reducing the amount of storage and the number 
of operations involved in a detection process of System 
B. On the receipt of r, the detector expands each of the m 
vectors (Q, .,] into the corresponding four vectors {P,}, 
where x, has the four values +1 and +3, and the 
detector evaluates c, (eqn. (12)) for each of the 4m 
vectors {P,}. The detector also expands each vector 
Q,- into four vectors {P,}, where x, has the four values 


The Radio and Electronic Engineer, Vol. 51, No. 9 


DETECTION PROCESSES FOR A 9600 bit/s MODEM 


tjand -3j( =, —1), and again evaluates c, for each 
vector P,. The detector then discards from each group of 
four vectors {P,} (originating from a single vector Q, . , 
via either of the two expansion processes) the vector P, 
with the largest cost c,, leaving three vectors in each 
group. Bearing in mind that yọ = I. this can be done very 
simply without in fact evaluating c, for any vector P,. 
but using just the real and imaginary parts of the 
quantity 


" n 
d, — Fi Y Nx s (n z X Xk 2 =x, (17) 
h-1 

in equation (12). for the different possible values of x,. 
Thus. from the values of d, for the six vectors {P,}. in the 
two groups of three vectors originating from any single 
vector Q, ,. the resulting values of d, for the nine 
vectors {P,}. given by all combinations of the three real 
and three imaginary values of x, in the six vectors. are 
very easily evaluated, to give the corresponding costs 
HULP} of the nine vectors. computed according to 
equation (13). x, is. of course, complex valued in each of 
these vectors. The detector now has 9m vectors [P,j 
together with the associated costs. The detector then 
determines s,. , from the value of x,_, in the vector P, 
with the smallest cost, and the detection process 
proceeds exactly as for System B. 


7 System E 

This is a simple modification of System D in which each 
of the m stored vectors {Q,_,} is expanded into four 
vectors (P,j. where x, has the four values +1 and +3. 
and also into four vectors {P,}, where x, has the four 
values +j and +3), exactly as before, but now the 
detector selects from each group of four vectors {P,} the 
two vectors with the smallest costs {c,}, the basic method 
of selection being that described for System D. From 
the four vectors [P,] derived from any single vector 
Q, ,. X, being real-valued in two of these and imaginary- 
valued in the other two, the detector forms the four 
vectors {P,} having the complex-valued [x,j given by all 
combinations of the real and imaginary values of x, in 
the original four vectors, and evaluates the associated 
costs {|U,|?}. The detector now has 4m vectors {P,} 
together with the associated costs, and proceeds with the 
detection of s,_,, and the selection of m vectors {Q,} 
exactly as for System B. 


8 System F 

This is a simple modification of System E which operates 
exactly as does System E except that the detector derives 
only three (rather than four) vectors {P,} with complex 
values for x,. from any single vector Q, |. ,. the vectors 
(P,! being those with the smallest costs (|U,[^] in the 
corresponding group of four vectors of System E. Again, 
the selection process for the {P,} can be implemented 
very simply. without in fact evaluating the costs. which 
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Fig. 3. Telephone-circuit 1. 


are determined after the selection process. Thus the 
detector generates 3m vectors (P,] together with the 
associated costs {/U,|?}, and then proceeds with the 
detection of s,_, and the selection of m vectors (Qj). 
exactly as for System B. 


9 Computer Simulation Tests 
The tolerance to additive white Gaussian noise of the six 
systems A-F and of a conventional non-linear (decision- 
feedback) equalizer have been determined by computer 
simulation over models of four different telephone 
circuits, using the arrangements shown in Figs. | and 2 
and described in Sections 1—8. In every case (including 
that of the equalizer) it is assumed that the adaptive 
linear filter has an appropriately large number of taps 
and is accurately adjusted to perform the ideal linear 
transformation described in Sections 1 and 2. 

Figures 3-6 show the attenuation and group-delay 
characteristics of the four different telephone circuits 
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Fig. 4. Telephone-circuit 2. 
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Fig. 6. Telephone-circuit 4. 


used in the tests, the telephone circuit forming a part of 
the transmission path in Fig. 1. Figure 7 shows the 
resultant attenuation and group-delay characteristics of 
the equipment filters, these being here considered as 
operating on the transmitted bandpass signal, rather 
than on the baseband signals in the transmitter and 
receiver as shown in Fig. |, to demonstrate more clearly 
the effects of the equipment filters in limiting some of the 
distortion introduced by the telephone circuit. The 
equipment filters are taken to include the filtering 
required to convert the sequence of impulses at the 
transmitter input (Fig. 1) into the corresponding 
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Fig. 7. Combination of transmitter and receiver filters. 


rectangular waveform used in an actual modem. 

Table ] shows the sampled impulse-response of the 
linear baseband channel, sampler and adaptive linear 
filter in Fig. 1, for each of the four different telephone 
circuits tested. This gives some idea of the resultant 
distortion introduced by each channel, bearing in mind 
that the ideal sampled impulse-response has the first 
component equal to unity and the remainder all zero. It 
can be seen from Figs. 3 and 7 that nearly all the signal 
distortion shown in Table 1 for the telephone circuit ] is 
in fact introduced by the equipment filters. The 
telephone circuits 1 and 2 introduce negligible and 


Table 1 


Sampled impulse-response of baseband channel and adaptive linear filter in 
Fig. 1, for each of the four telephone circuits. 


TELEPHONE CIRCUIT TELEPHONE CIRCUIT 


H 2 


REAL IMAGINARY 
PART PART 
1 
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PART PART 


TELEPHONE CIRCUIT TELEPHONE CIRCUIT 


3 4 
REAL [IMAGINARY REAL [IMAGINARY 
PART PART PART PART 
1.0000 0.0000 . 0000 0.0000 
0.4861 1.0988 | 0.2544 .9941 
-0.5980 0.0703 -1.7394 -0.2019 
0.1702 -0.1938 | 0.6795 -0.8086 
-0.0245 0.1000 | 0.0408 0.5113 
0.0100 | -0.0258 | -0.1189 -0.1463 
-0.0134 0.0110 0.0343 0.0420 
0.0056 -0.0042 -0.0185 -0.0364 
0.0003 .0003 0.0139 0.0216 
-0.0008 0041 -0.0102 -0.0009 | 
. 0000 | -0.0061 | -0.0019 -0.0034 | 
.0007 -0.0007 0.0037 -0.0046 
0.0037 | 0.0002 -0.0028 -0.0006 
-0.0019 | -0.0025 | -0.0019 -0.0046 
0.0020 | 0.0008 | 0.0083 .0022 
0.0005 | -0.0002 |-0.0056 0059 | 
-0.0022 0.0002 | -0.0046 -0.0028 
0.0007 | -0.0005 0.0049 -0.0019 
-0.0008 | 0.0002 | -0.0009 0.0037 
0.0005 0.0005 -0.0009 0.0003 
0.0000 0.0000 0000 0.0000 
0.0000 oO. . 0000 0.0000 
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typical levels of distortion, respectively, whereas the 
circuits 3 and 4 are close to the typical worst circuits 
normally considered for the transmission of data at 9600 
and 600-1200 bits/second, respectively. The telephone 
circuit 3, which is close to the Post Office network N6, 
introduces severe group-delay distortion as well as 
considerable attenuation distortion, and the telephone 
circuit 4, which is close to the Post Office network N3, 
introduces very severe attenuation distortion. 

Figures 8(a)-(d) show the performances of Systems 
A-F and of a conventional non-linear equalizer, for the 
telephone circuits 1-4. The signal/noise ratio is taken to 
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be y dB, where 
Vy = l0log,o (E/4No) (18) 


E is here the average transmitted energy per data-symbol 
s;, at the input to the transmission path (Fig. 1), and No 
is the two-sided power spectral density of the additive 
white Gaussian noise at the input to the receiver filter. 
The 95% confidence limits of the curves in Fig. 8 are 
generally better than +4 dB. System A has been tested 
with both m = 16 and 32 and with n = 7, where m is the 
number of stored vectors {Q,! and n is the number of 
components in Q,. Systems B-F have all been tested with 
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(d) Operating over telephone-circuit 4. 


Fig. 8. Variation of error rate with signal/noise ratio for the different systems. 
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m= 6 and with both n = 3 and 7. For the telephone 
circuits 3 and 4, System B has in addition been tested 
with both m = 4 and 8 and with n = 7. 

The phase of the sampling instants {iT } at the receiver 
has a much more significant effect on the tolerance to 
noise of the conventional equalizer than on that of any of 
the other systems tested. The sampling phase used over 
each channel has therefore been selected to give a good 
performance of the conventional equalizer, with a 
tolerance to noise within about } dB of optimum, which 
is probably as good as would be achieved in practice. 

It can be seen from Fig. 8 that the Systems A-F all 
have a much better overall tolerance to additive white 
Gaussian noise than the conventional non-linear 
equalizer. The System À with m — 32 (32 stored vectors) 
and n = 7 (a delay in detection of 7 sampling intervals) 
has the best performance of all systems tested, as would 
perhaps be expected, but when m= 16 and n = 7 its 
tolerance to noise at the lower error rates in the {si} 
becomes inferior to that of the Systems B-F, with m = 6 
and n — 7. The latter arrangements all have a similar 
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Fig. 9. Variation oferror rate with an innaccuracy in the estimate made 

by the receiver of (a) the level and (b) the carrier phase of the received 

signal, for the System E with m — 6 and n — 7 operating over the 
telephone-circuits 1 -4. 


tolerance to noise at the lower error rates, this being 
significantly better than that with m = 6 and n = 3, over 
the telephone circuits 3 and 4. No very useful advantage 
in the tolerance to noise of the Systems B--F seems likely 
to be achieved by increasing m to a value greater than 6, 
but when m is reduced from 6 to 4 there is a noticeable 
degradation in performance. At error rates around | in 
10°, the Systems B-F, with m = 6 and n = 7, have an 
advantage in tolerance to additive white Gaussian noise, 
over the non-linear equalizer, of approximately 1, 14, 3 
and 51 dB, for the telephone circuits 1, 2, 3 and 4, 
respectively. 

Figure 9 shows the effects on the error rate in the {sj} 
of inaccuracies in the estimates made by the receiver of 
the level and carrier phase of the received signal, when 
System E with m 2 6 and n=7 operates over the 
telephone circuits 1-4. It can be seen that, over the 
telephone circuit 3, the system can tolerate an inaccuracy 
of about 1 dB in the estimate of the received signal level 
or else an inaccuracy of about 2 degrees in the estimate of 
the received signal carrier phase, for an increase of two 
times in the error rate of the {s;}, this representing a 
reduction of less than 4 dB in tolerance to additive white 
Gaussian noise. It seems likely therefore that in a well- 
designed modem and with the correction of the more 
serious phase jitter in the received signal carrier, no very 
serious degradation in tolerance to additive noise should 
be experienced due to the inaccurate estimation of the 
received signal level or carrier phase. 


10 Conclusions 

Of the various arrangements tested. Systems E and F 
with m — 6 and n 2 7 appear to be the most cost- 
effective. They promise to give a substantial 
improvement in tolerance to additive noise over the 
conventional non-linear equalizer, without involving a 
great increase in equipment complexity. Bearing in mind 
that the telephone circuit 4 is basically similar to the 
typical worst telephone circuit likely to be experienced 
over the switched telephone network in this country, it 
seems that so long as the accurate convergence of the 
adaptive linear filter can be achieved and appropriate 
steps are taken in the design of the modem to combat the 
effects of carrier-phase jitter, correct operation of the 
system should be obtained over nearly all circuits on the 
public switched telephone network. 
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